PICS e eR ele 


W'9 PHL G HT SRM CHE ETS DS eT ERE EYE 


For Reference a a ie 


NOT TO BE TAKEN FROM THIS ROOM 


eke ures Cee ti 


(aevlotsver 


Gx ipnis 
UNIVERSITAT 
LABERTAEDSIS 


on 


haw 


“= ‘ : ane, = oe os 9 

ides am Bic at 
he = : a x VF ‘ei \ 

y } i taille A _ we ¥ } 
Se aa 


i a 7 A ae ; . a Viva 
4 


A w 4 sa ie tt er ay Cr 
SEY af eis +d a 2 £ 


: irae 
' 7 wa 


: sR 


: f 
Zh _ 2 tee bo ee 
et ES : ' af ee a 


rer AZ) at 


¢ Ye Y, 
= 
‘ 
a ( A ie «+ 7 
$2 ' if 2 here 
' ee: 
: j : 


ee es Lie we a 
~ | ne - ‘ah 
oe" si ‘ ou) c 2 eee Ar si | 
| ear 
Seimei eéeved BBC BO! ru iP 
ae a i, 
Tt: Ge oc Ys, ear | a ies 
/ ' ae es ah \ ne 
jey¥oc bite Sewhs) Gabe Ties so, Peet ve sbia Ne a 
i 2 at ic) a as 
sui ¥Taeeaken i S)aeaki« i. atearal tar 
a | os ae 
To ete waa * ay | ene 
v. Metees CRAG ease perieaip niieus Oa ee , 
hee Be ¢ sk 4s) . a 
a * DI ry ve a6 me fi Say Fs Pr “male A ls 


‘ AG -y 
patwe’ SAS Supiisis ns: a Ai aq ra Stent ds. oT er , 
“ 7 saa canal natah 


ee 


: aa 1) iB: Me 
y > Si5 am, - b eitiy i) 7 - 41% : 4 


daaea eis 
“atte a oo 


; Seyemets, Seok: 


buat ta wey, 


i aa 


THE UNIVERSITY OF ALBERTA 


DHEPDIEUTE XY MODEL 


by 


Ce) Correia LEP Rene REE ye 


A THESIS 
SUBMITTED T0 THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
Phe PARMA) FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE ‘OF 
DOCTOR OF PHILOSOPRUY 


DEPARTMENT OF PHYSICS 
EDMONTON, ALBERTA 


FALE, 1976 


Digitized by the Internet Archive 
in 2023 with funding from — 
University of Alberta Library 


https://archive.org/details/Reeve19/6 


To my wife 


Monica 


iv 


ABSTRACT 


The critical properties of the quenched and annealed 
site diluted spin-5 XY models are investigated, for three 
cubic lattices, using the finite cluster method. .The lines 
of second order transition points for both models in the 
temperature-density plane are located. Benavior of the 
€ritical exponent y is: found not to. agree with the ideas of 
universality. The annealed site XY model is shown to be 


equivalent to the Takagi model of He>-He@ 


mixtures. °‘Contours 
of constant reduced density n = na + Ny in the TX, plane are 
determined and whe EVICEeNa Cam POINtLLS enproximacciwy: |ocated 
for the f.c.c. lattice. ._ The two dimensional quenched and 
annealed site XY models are analysed for the triangular 


lattice. The resuits seem to indicate that a phase transition 


does occur. 
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CHAPTER I 
INTRODUCTION 


Many Ait serene typeswvot, pnasectransiitam occurain 
Doture. put 1¢ 1S onlyvin recent years that thereihas been 
a concerted effort to understand these phenomena collectively. 

Van der Waals formulated in 1873 the jaw of 
corresponding states, which asserts that the equation of 
State for all fluids has the same form. This theory has 
enjoyed considerable success and it has served to describe 
a number of liquid gas transitions reasonably weil. 

Magnetic systems however, could only be explained 
from a quantum mechanical viewpoint. In simplistic terms 
the Pauli exclusion principle tends to keep electrons with 
parallei spins apart and in doing so reduces their Coulomb 
repulsion. The difference in energy between the paraliel 
and antiparallel configurations is the exchange energy. If 
the increase in kinetic energy associated with a parallel 
alignment is less than the decrease in potential energy, 
then magnetic ordering can accur. The physical mechanisms 
involved in magnetism are discussed in detail in Mattis 
(1965). 

One of the simplest, yet most successful moaels 
of magnetism is the Ising model. The free energy of the 


Ising ferromagnet on the square lattice was found 
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analytically by Onsager (1944). As yet no exact results 
concerning the phase transition for the three dimensional 
Ising model are known. Because the Ising Hamiltonian in 
the Pauli representation is diagonal, it is often referred 
to as a classical model. The XY model, which is a special 
case of the anisotropic Heisenberg model is “fully" quantum 
mechanical. The XY model has only been solved exactly in 
one dimension. 

Since exact solutions are out of the question for 
most realistic models of magnetism a number of approximate 
metnods have been developed. The principle closed form 
approximations are the metnods by Weiss (1907), Bethe (1935), 
Bragg and Williams (1934) and Kikuchi (1951). Also of 
importance are the Green's functions techniques. These have 
been reviewed by Zubarev (1960) and Tyablikov (1967). 

By far the most successful method of elucidating 
the critical properties of direct exchange type modeis is 
the method of series expansions. In this method, thermo- 
dynamic functions are expanded about the ordered state or 
about the disordered state in powers of a suitable variable 
such as temperature or density. The series are then analysed 
using various methods to reveal any singular behavior which 
may indicate the occurrence of a continuous (or second order) 


phase»transition. It. is this method we exploit here. 
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Before briefly describing the purpose of this 
thesis though, we mention the very recent development of 
renormalisation group theory. This is an approximate method, 
but is enjoying considerable success as an approximation 
applicable to classical models. For reviews the reader is 
referred to Ma (1973), Wilson and Kogut (1974) and Fisher 
(1974). In its original form renormalisation group calcula- 
tions are done in a momentum representation, but Niemeyer 
and Van Leeuwen (1974) have developed the theory in 
configuration space. This newer formulation appears to 
be conducive to calculations on quantum systems (Rogiers and 
Dekeyser (1976) and Betts and Plischke (1976)). 

Comprehensive reviews covering most aspects of 
critical phenomena are contained in the set of volumes 
edited by Domb and Green (1971) and the articles by Domb 
(1960) and Fisher (1967). De Jongh and Miedema (1974) have 
extensively reviewed the experimental situation. 

The XY model can be explained from a physical 
standpoint as follows. Consider the following model of a 
magnetic crystal. At each lattice site is located an 
elementary magnet which is allowed to rotate in a pre- 
determined plane. These magnetic moments can be represented 
by two dimensional vectors. If the interaction between 
neighbouring sites is proportional to scalar product of the 


vectors representing the magnetic moments at those sites, 
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then the model is called the plane rotator model. Further, 
if the elementary magnets are allowed to move in three 
dimensions but are still restricted Lorine naee in two 
dimensions, then the model becomes the classical or spin 
infinity XY model. 

The spin-y XY model is the quantum analogue of 
ere Spin intinity  X¥ "model ?“"Theldifferenceis that the 
elementary magnets at each site are spins which interact in 
two dimensions only, while there are still three rotational 
degrees of freedom available to each spin. It is perhaps 
more enlightening to regard the spin-5 XY model as describing 
a Heisenberg magnet in which the crystal field has the 
appropriate symmetry and strength so as to prevent inter- 
action between spins in one particular direction in spin 
Space. Several examples of the spin-z XY model have been 
identified experimentally and are listed in Chapter III. 

| The’ eraticatsproperties of*the spin-z XY model 
have been considerably elucidated by Betts (1974) and co- 
workers and by Dekeyser and Rogiers (1975). 

OnevoTt’the systems studved"inw thts’ thesTs 1s" the 
quenched site spin-y KYSnOdete, ints’ ts asModeleror a=: dir tye, 
crystalline magnet in which the impurity atoms are frozen, 
in random distribution, into the lattice. Study of such 


a system is motivated both from a theoretical and experimen- 


tal viewpoint. 
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Theoretically there has been considerable specula- 
tion as to the nature of the phase transition of diluted 
magnets. It has been suggested that phase transitions in 
random systems may not occur at a well defined temperature 
but over a small finite range of temperatures. Also there is 
the question whether or not scaling theory and the univer- 
Sality principle hold for such systems. These questions are 
expanded on in Chapters III and V, where extensive references 
can be found. 

To date experiments on random systems have been 
restricted to mixtures of two different types of anti- 
ferromagnets. However de Jongh (private communication) has 
begun studies of some randomly diluted XY antiferromagnets 
and anticipates that he will be able to determine the criti- 
cal temperature for a number of different impurity 
concentrations. 

Besides being a model for a crystalline magnet the 
XY model can also be used as a lattice fluid model for the 
normal fluid to superfluid transition in He’. Consider the 
He? fluid as a.collection of bosons arranged randomly in 
space. If we then break up the space into cells large enough 
only to be occupied by a single particle, then some of these 
cells will be occupied by Hes atoms and some will be vacant. 
The XY model represents such a collection of bosons and 


allows the Het atoms to move to unoccupied neighbouring cells 
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and in this way accounts for the kinetic energy of the system. 
Theoretical results for the critical properties of the XY 
model. have been compared with experimental data for the 
normal fluid to superfluid transition in Hes with excellent 
agreement (see Betts (1974)). The lattice fluid picture of 
the XY model is described in much more detail in Chapter III. 
This brings us to the second model studied in this 
thesis. This is the annealed site spin-5 XY model which 
corresponds to He@ with classical impurities added: ‘Such ‘a 


model could conceivably be used as a model of He 
mixtures with the kinetic energy and fermi statistics of 
the He? atoms being disregarded. This approximation is 
presumably best for temperatures above the tricritical 
temperature where phase separation first occurs. The 


annealed site XY model of He® ~ He’ mixtures ignores all 


oer bees He? - Ta and He* - Het 


but the hard core He 
interactions. However the critical properties of interest 
to us arise from the kinetic mobility of the He’ atoms and 
neglect of the potential tail will not alter the values of 
the critical exponents. This invarience of critical proper- 
ties with respect to perturbative effects arises from the 
universality principle as stated in Chapter II. 

The aim then is to find the phase diagram for 
the spin-+ XY model in the temperature-density plane and 


2 
compare this with the experimental result. 
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The plan of the remainder of this thesis is as 
follows. In Chapter II some definitions and phenomenalogical 
theories of critical phenomena are outlined. The XY model 
of both magnets and lattice fluids is introduced in Chapter 
III. Also in Chapter III we define precisely the models 
studied, namely the quenched and annealed site spin-z XY 
models. The finite cluster method of constructing series 
expansions is described in Chapter IV. In Chapter V the 
results of the analysis for the quenched site model in three 
dimensions are presented, while the results for the annealed 
site model in three dimensions are given in Chapter VI. The 
two dimensional quenched and annealed site models are 
discussed in Chapter VII.- Final conclusions are given in 


Chapter, VITI-. 
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CHAPTER II 
CRITICAL PHENOMENA 


Sel aelavLRoduction 

The definitions and general results from the 
theory of critical pnenomena that are referred to in later 
chapters are outlined in the following sections. 

Basic definitions of some thermodynamic functions 
and their critical amplitudes and exponents are given in 
PpecturOn c.i. Ie tneorye oT itricritcal scaling is briefly 
described in Section 3... Section 4 deals with the Fisher 
renormalisation of exponents. The concept of universality 
of critical phenomena is discussed in Section 5. 

Throughout this thesis we will use mainly the 
thermodynamic notation appropriate for magnetic systems. 
Tne conversion from magnet to fluid language for any 
formula is usually easily facilitated by making the 


substitutions 


te wie 


(2esl231)) 


phe kat 


where M and H are the magnetisation and tne ordering 


magnetic field respectively. The density and chemical 
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potential of the fluid system are labelled p and u 
respectively and have tne values Po and He at the critical 
point. When reference to phase separation in fluid mixtures 
is desired, tne concentration of one of the species is 
related to tne magnetisation and tne difference in the 
chemical potentials of the two species is the ordering 


he Bee mixtures the order 


r1erd.. in particular for) He 
parameter is the He® concentration X 3 and tne ordering 
field is U3 - Ug» the chemical potential difference. 

We leave until the next chapter a discussion of 


just now the dilute XY model can be used to test the 


phenomenological and thermodynamic theories given below. 


eee  Cribrcalwexponents 


Tne thermodynamic functions usually dealt with 
when investigating a phase transition are the order 
parameter and various response and correlation functions. 
bet F{(T;sH) be the free energy of a system, where H is the 
ordering field and irrelevant variables have been 


Suppressed. Then tne following functions can be defined: 


the magnetisation 


M(T,H) = (9F(T,H)/3H)> (aso 1%) 


the isothermal susceptibility 


x(T,H) = (abi(T,H)/9H), (C228) 
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and tne specific feat 


COURS MER SSE (THY STS) 


‘ (2.9233) 


Higner derivatives of the free energy are frequently 
encountered in tne tneory of critical phenomena and are 


usually denoted 


FS) (TH) = (aE (TH) /aH*) (2.2.4) 


T 


In addition the spin-spin correlation functions 
characterise the microscopic behavior of the system. The 
pair correlation function of greatest interest for the XY 


model is 
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wnere in (2.2.5) a isethe x component of Unerspin at 


@elatttce site a’ distance r' from>some reference lattice 
Site-ae'o. The*quantity S ts the magnitude ofthe spin dt 
each s#te. For the XY model <S*> is the order parameter, 


SO°FOrRT. > li (the critical temperature) the self correlation 
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tions the reduced susceptibility may be expanded as 


term <S is ‘zero. In terms of the pair correlation func- 


Pens Loh) (2.42.6) 
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In general the reduced form of any thermodynamic quantity is 
formed by dividing that function by the same function for a 


paramagnet at the same temperature. Thus we have 
m(T,h) = M(H,T)/M(H+@,T) i aeeerya) 
for the reduced magnetisation and 


C = €,(T,H)/Nk T G27) 


H B 


is the reduced specific heat. The reduced susceptibility is 
i Miles ewe: C1, Hy VT SH) ee 


where h = HA/KeT is the reduced ordering field and u the 
magnetisation per particle. 

The modern theory of critical phenomena is 
formulated in terms of critical exponents and critical 
amplitude ratios. Let t = Or ae be a measure of the 
distance in temperature from tne critical point. Tnen the 
Criticw#l-exponent for T > ips and h = O for any thermodynamic 
quantity f(t) is 
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implied asymptotic form for f(t) is 
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where A is cailed the high temperature amplitude of f and 


(fot OYE oil. J.) . nes usial notation for, the: moremcommon 


: 
exponents and amplitudes is as follows: 


m(t,0)  B(-t)® t<0 
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Get at la! t<0 
acts moh eee 
ct t>0 


and P 
“Ay 


uel ee0) athe) (2) ashe ee a (tao) 


(2) ene ooea( 22) 


F(22) 4 OV) DE £ Gte.0}) 


The correlation function has the form 
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and the correlation lengta k(t) has the asymptotic form 
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Parean extensive list of critical exponents and‘ their 
values for some experimental situations the reader is 
pererred to Stenhey (1971). For the pure spin-z XY model 
the results of Betts (1974) and co-workers and Dekeyser 
and Rogiers (1975) indicate that aw0 (a near logrithmic 


divergence) and y = 1.333 + 0.002. Other exponents may be 


obtained from scaling theory as explained in the next section. 


Boe oCaling,.neory 
ee or bencalepoint scaling 

Scaling theory is a phenomenological theory of 
the behavior of thermodynamic functions near a critical 
point. The theory was originally formulated by Widom 
(1965a,b) and Domb and Hunter (1965). Kadanoff (1966) nas 
given heuristic arguments in support of the theory. An 
outline of critical point scaling is given in Stanley (1971) 
and extensive reviews by Vicentini-Missoni (1972) and 
Levelt Singers (1974) record much of the known experimental 
evidence in support of scaling. 

The, static scaling hypothesis -assertsstnat tthe 
singular part of the reduced free energy near the critical 


region is a generalised homogeneous function. This is 
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usually expressed matnematically as 
Pots ih = en een) (223hap 


where the exponent A is the gap exponent whicn is equal 
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By taking various derivatives of the free energy 
equation (2.3.1), every exponent previously defined may be 
written in terms of the two exponents A and 8. As a 
consequence of scaling, Rusnbrooke's (1965) thermodynamic 


inequality 

Ge ite2eet 7 oie ee C2e352) 
and Griffith's (1965a,bd) thermodynamic inequality 

BcOpert ht (28264 Rees 


are both satisfied as equalities. Scaling theory also 


predicts that corresponding high and low temperature 


4 


exponents are equal, for example y = yy Che aek 
Experimentally the scaling hypothesis can be verified 

for many substances and for several different types of 
second order transitions. In particular the: scaled 
magnetisation, m(t,h)/(-t)®, has been plotted against the 
scaled magnetic field, WE aac by Kouvel and Comley (1968) 


for nickel (metallic ferromagnet) and by 40 and Lister 


ps 7A ‘ aT, / 7 ’ ui 7 ru aah - We tA - 
4 et iin Pm | 4.6 — ars : a ic i i 


o on a | aa 


“ ‘ : a MP as. i bp aly < 
= .) in 
aes t ay, 4 | +p ro S¥ nj iN erg AY id 2ub net) ok 
ae We : rt, a D ae aun Be es 
eT Sere win by D1 Dh) Pere of raney a a “yi a é . é , 


he hie! e ¢ paved: a) x: bey ni 
f | ya} ; kad {aii hiss 


Prana bors sor Cent ‘4 ADO7 20 Sit EH ths e | 


} se 
2 | n 
ae hs toy 
~ 
, ; 7 Ga t a 2 
feo Po.’ 
a ‘ . + 


vegies phe os aM 257 yet tee 
3 ’ * el ai ; ~ ; 
Soinrs3 oa i= i) Te» th As an gree bah. +2) ase 


; | a 
ae: Sie a a Ae SUuapanen J tige ALINE. 
Beit itav 20° We zi railad ; nit ng) S| 
Pal ne ier ; 


i? RRNA S RAialy 1p | seas “oA Way 
a patane ane i, ities a 6; batats i 
he (ag daeitgelbeisaty nbeige Be R "Ga Nig 
ngs) eon 7 tedugi knees ie ; 
P og Taed bi eyugs, Sar PAs, ha Pa 


5 ‘ en pakeny ty a - ie dorsi 1 ban 


5 


(1969) for CrBr., (insulating ferromagnet). Their plots show 
ehat fat: thie pOINLS peor 1) < Li iO mOn CurCU RY ema tidisd || 
Doarmts Tor: [> ihe fit on anotner, verifying the predicted 
temperature independence. 

As stated above only two exponents are needed to 
obtain the complete set of exponents for thermodynamic 
Guanuwities: For the XY spin-z model, assuming y=4/3 and 


a=0, scaling theory predicts y=2/3, A=5/3, 8=1/3 and n=0. 


Zooey Un icriticalescalings theory 

Petrie yh ised limp Onn ty waSeGetTdineduMb year t t1 tis 
(1970), is characterised by the existence of two competing 
order parameters which simultaneously become critical. 
Examples of tricritical behavior can be found in metamagnets 


like FeCl, (Yelon and Birgeneau (1972)), Ni(NO,),°2H,0 


2 
(Schmidt and Friedberg (1970)) and dysprosium aluminium 
garnet (Landau and Keen (1972)). Some systems which display 
structural phase transitions such as NH, Cl (Garland and 
Weiner (1971)) aver tricriticalapoints. Of particular 
inteyest tous i1s.the tricritical’ behavior. of tes elhes 
mixtures because, as we show later, the annealed site XY 
model is conceivably a good model for this system. 

Figure 2.1 shows the phase diagram for He? - He’ 
mixtures in the TA plane, where A is the chemical potential 


difference Uz - Hy 
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Figure 2.1 Temperature T versus chemical potential 

difference A for He*?-He* mixtures (qualitative). 
Tne heavy solid curve represents the second 
Order phase transition line wnile the dashed 
curve represents the first order transition 
line. The scaling axes h and g are plotted and 
the curve g=h® witn ¢ > 1 divides the plane 

into regions A and A™ (u-like exponents) and 

and B“ (t-like).exponents. 
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morrconstant A :< A, tne transition is from the Hien ae 
fOmcne normal fluid, state. of He’, When A > A, the 
transition is first order,and the line defines the phase 
boundary of the He'-rich and tne Hee-rich phases. 

Recently the homogeneity hypothesis has been 
exrendeducouinciludeytricritical. points. Riedel (197.2. 1974), 
Griffiths (1973) and Hankey et al. (1972) nave developed 
Eneatneory of tricritical scaling. Let, be the: superfluid 
Orderiparameter, and, < the field conjugate to w. «Tne field 
PecuyS UCd hive lnrealisabile.. ‘Figure 2,1, 1s the Hee-He@ 
phase diagram in the c=0 plane. In the field space 
(jAvplane) the: tangent to’ the line.of critical points 
Singles out a special direction whicn is used as one of 
the scaling directions. The co-ordinate along this 
direction is denoted as h. Anotiner co-ordinate; g, can be 
chosen so that it makes an arbitrary non-zero angle with 
Prev iiedi Gecticon. sled loUure 2. lonvand gare. drawn at 
Pigntwangiess. one third .scaling.variaole Isic. 

Tne tricritical point homogeneity hypothesis can 


tnen pe written as 


AU Ue a a fonegece (2534)) 


waere i i.G,c)vis again, the most. singular part of the 


reduced free energy. Griffiths (1973) pointed out that the 
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above assumption results in two sets of tricritical exponents 
depending on path of approach toward the tricritical point. 
The two sets are distinguished by subscripting the exponents 
mirieaet On-u. In (2.3.4) the t.subscripted exponents are 
used and » is called the cross-over exponent, the value of 
which determines tne different scaling regions. Throughout 
the remainder of this discussion we assume that » > 1. For 


He--He? 


mixtures experimental evidence (Godner et al. (1974)) 
indicates that @~ 2. 

Atytnts *goint' we make “a stignt*digression to 
describe some results obtained by Riedel (1974) and co- 
workers using renormalisation group techniques on tne 
Blume-Emery-Griffiths (BEG) (1971) model. Tne BEG model 
is essentially a lattice fluid spin-1 Ising model of 
He?-He- muxcuresr “kn this! model. a He” barticie,. oOnsene 
ith lattice site, is represented by the Se =e S teat Gieaed 
He? particle by tne S<=0 State and the absence of a 
particle by the Se = -] state. Tne value o=2 was obtained 
by Riedel and Wegner (1972) for the BEG model. Other 
exponent values they obtained will be presented later 
Tn thise@ecti on. 


The role of tne cross-over exponent is illustrated 


below. by example. Assuming the ansatz (2.3.4) and using 
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Putting c=0 and X = h 16 G2 345 YaGives 


—o(2-a,-A 


v(h,g) = h t) 


v(h/|h] ,g/n?) oaean) 


Equation (2.3.6) iS-walid when g « h? which we define as 
the region for the u exponents and oe o(2-a,-A,). Now 
putting c=0 and jeg i/o Tua 2 33.5) Ouves 
-a,t+2-A 
Drege! Sy(n/g!/®,g/| 9) (2.3.0 ) 


Pauation) (223.7) issvalid when g > n? which is the region 


for the t exponents and 


By By ea- apt A, z B/¢ (2753 58)) 


In general the t exponents are related to tne u exponents 
by a factor $9, except for the specific heat exponents where 


the relation is 


a, 7 oe o(2-a, ) ero) 


The path g=h? is sealed!’ thenscalingeboundarye) 71h figure 
2.1° thes line g=h? issplotted for deetd andathe different 
scaling regions are identified. 

Quantities expected to be singular at the tri- 


critical point only are the concentration X 3 and the 
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concentration Susceptibility defined by 


W = (8X3/A)_ 


where X, = (3F/9A), 
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If X, is the value of X, at the tricritical point, then 


using the now standard notation we have 


W 
is C ! g>h? 
ine en OA! GD 
Sw a3 + u g<h? 
and -i 
c . g>h? 
Ww { -d 
be g<h? 


An analogous procedure to that given above gives 


a r a o(1-a, ) 


and r A, = oa, 
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Experimental results for Heo-Hes mixtures by 


Goellner et al. (1973) and renormalisation group calcula- 


tions on the BEG model (Wegner and Riedel (1973)) yield the 


classkeal set of tricritical exponents, namely ¢=2, a,=1/2, 
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Next we look briefly at cross-over effects. 
Cross-over arises from the possibility that there exists 
mregronewnere both critical and tricritical scaling Jaws 
hold. 

Consider the susceptibility in the region near 
both the tricritical point and the line of second order 


piace transitions (irraicritdcaly scaling predicts, that 
oY, $ 
VEO Nee AL XAG SAH) (253014) 


and critical scaling asserts that 
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The result is that the amplitude of ee varies as 


ON tee) 


Pee an (2773) 86) 


which diverges or converges as h moves closer to zero, 

depending on the sign of Oe In the cross-over region 
the observed exponents (say y for example) will not have 
ehinerecnitical or tricritical values, but some effective 


value 


which depends on g (Riedel and Wegner (1974)). Wegner and 
Riedel (1973) nave also found that at the tricritical point 
there may be a logrithmic correction to the simple power 


law behavior for some thermodynamic functions. 


2.4 Fisher Renormalisation 

This section deals with the renormalisation of 
exponents by "hidden" variables. 

The results we quote. below are for thermally 
diluted systems, where the impurity concentration is the 
"hidden" variable. The renormalisation formulae were 
derived in general by Fisher (1968) and are summarised by 


the formulae 
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where ey is the exponent y when the system contains a 
concentration x of thermal impurities. All other exponents 
are renormalised as y and 8. The pure exponent values are 
found only when x=0. The formulae (2.4.1) are subject to 
tre cond? Cionmthtarbs: ot 10.2 Iaeethe event: that) ai Oi thiat 
is the specific heat is discontinuous, but finite (or 
possibly non-analytic) at ae then all exponents vis BY. 
etc. will be observed to be greater than their values of 
x=0. Only when the specific heat has a logrithmic diver- 
gence at ay is there no change in the observed value of 
exponents from their "pure" values. 

The formulae (2.4.1) have been rigorously shown 
to hold (Essam and Garelich (1967)) for the Syozi (1965) 
model, which is Beenie to an annealed bond spin-y 
losis ysite mea lin the: caserof they modell cine:rspeicitfirc 
heat singularity is logrithmic or nearly so. This means 


that the Fisher renormalisation of exponents will be 


undetectable even if it is present. 
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Pepa wuni versal ity 

The universality hypothesis has been stated by 
Kadanoff (1971) as follows. 

"All phase transition problems can be divided 
into a small number of different classes depending on the 
dimensionality of the system and the symmetries of the 
Ordered state. Witnin each class, all phase transitions 
have identical behavior in the critical region, only the 
names of the variables are changed." 

The concept of universality is inspired in part 
by experimental evidence. Most simple fluids have a 
Similar set of exponents for the liquid-gas transition 
which are very close to the 3-D nearest neighbour spin-s 
Pree xponenitwul eC. yo he 25 ‘ae 0. h25 ‘and @ 2 OU3a25 
Many different magnets also have similar exponents. 
Specific examples can be found in Kadanoff (1971) and 
Stanley (1971). 

In addition to the dimensionality of the system 
and the symmetry of tne ordered state the universality 
classes are thought to depend on the range of the interaction 
and possibly for magnetic systems the magnitude of the spin 
per lattice: site... One of the’ earlieste formujations ofthe 
universality principle is due to Betts, Guttmann and Joyce 
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Singular parts of their free energies near their respective 


transitions satisfy 


=nf 
ny ybtyehy) om 57 ae) 


with Geis Gt, and ae = rats 2 ae) 


Where t = Tele U5 h = mH/kT and oF and Acie Crit cay 


scale factors. An immediate consequence of (2.5.1) and 
tee og eyensethat the critical exponents of) X and Y are 
equal. For the spin-spin correlation function the 


universality relation is (Ferer and Wortis (1972)) 


(Rta her) (265.3) 


ONES mah Bacal 


with R/O yey = RY /5 hy vasa se) 


where RY is the relative position of the two spins, jy is 
the nearest-neighbour separation and Ly 1s? arthird seattle 
Factor: 

Recently the’ concept of universality has been 
refined!by two-additional postulates. Firstly Stauffer, 
Ferer and Wortis (1972) and Ferer and Wortis (1972) assert 
thatetheimost singullar*pant of the free. energy’ int a*valume 


with diameter equal to the coherence lengtn is universal. 


As a consequence of this condition only two independent scale 
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factors are needed. Secondly Betts and Ritchie (1975) 
suggest that the most singular part of the energy per bond 
length is universal. This conjecture necessitates only 
one scale factor. However the one scale factor universality 
does not hold for the spherical model of Berlin and Kac 
(1952) in the same form as for other models. 

Some theoretical aspects of universality have 
been reviewed by Stanley (1974). In general the hypothesis 
as given above is consistent with a wide range of experimen- 
tal and theoretical data. However there is considerable 
doubt as to the applicability to systems having more than 
two spin interactions. The best example of this is the 
eight-vertex model, solved exactly by Baxter (1971, 1972). 
The Hamiltonian for this model contains four-spin interac- 
tions and the critical exponents depend on the coupling 
strength. This phenomenon appears to be a violation of 
universality, but Kadanoff and Wegner (1971) have shown 
that the order parameter does have some peculiar symmetries 
which indicates that the eight-vertex model is in an 
unusual universality class. However, Suzuki (1974) has 
proposed an alternative remedy by using a "weak" universality 
principle. Weak universality requires that all exponents 
should be defined in terms of the correlation length k and 


Nnoftat. she result 1s. that the’ ratio for any two exponents 


Lent “4 
» 
4 ra 
‘ yee | ; 
> Ar 
Le 
ss : ae “A 
” * é - * ¥ mri, & a . i wis 
‘e i re rik | 
d 
Z : a i! | 
j 
1 7 wy Sia tas 
4 c ph ‘ | i) 
~ { 1 -? ol By 
: at 
i, 
- sf " 
; ' 
‘ < 
~ 
, 
a7 Vor oe 
| .. #; 
- Meta 
7 ‘ 
f y=: 
a Bs 
yh ale. 
i ih SV 
. 
OF US OP hans tie 


: SNS Beast 
ry : q 
« ' pe t ‘ 
¢ Bach \ hohe ye 2 OF. ty 
= it > ie 
ce i ” ae | ; ; Pie 
‘iy & > on AO 4S Lar) ee 
; “9 a 10 — x 
“f : a 
Bi Gl ee tas ph agRe ats OHO cd 
. aa! | OL ae 
ie sei ae tee 
yen OR Olt ase tine Ansaptnaily 


St. BOD BVEe sunt 38 mea 


; ': J ‘a Sy? j i - 
wesw” oh bi Pails Veet an avd ten a, 
eee SI abe al 4: 


Cin 


for systems in the same universality class are equal. More 
particularly the values 6=5 and n=0 are universal values 
for three dimensional systems and y/v, B/v~, etc. are 


invarient within each universality class. 


CHAPTER iI! 
THE XY MODEL AND MAGNETIC DILUTION 


See A introductacon 


Before motivating a study of the XY model let 


lisautians ti <define what +t as. 


The anisotropic Heisenberg model for spin-+ 


2 
systems can be written as 
2,-1 8 eps yh OY II Zaz 
H = -(4S a EER, ES URASIO Ne a BOCAS AR ie fad Ea Gs SI: 
ae ey Oh ea 
-] Xx Zee 
-S : tm, Si + mH Sa (3m lee li) 


the indicesi andj run over all lattice. sites and Ss; 
is the a Cartesian component of the spin vector S of 
magnitude s. The exchange "integrals" re 
measures of the spin-spin coupling strength and the 


II 
and Jj are 


external fields H* and H% are coupled to the spin system 
via the coupling constants my and m, 
Three special cases of (3.1.1) can be identified. 
These are the Ising Hamiltonian (when jt = 0), the 
isotropic Heisenberg Hamiltonian (when gl = g+) and the XY 
Hamiltonian (when gil = 0). Both the Heisenberg and Ising 


models are well established as models for certain strong 


ferromagnets and antiferromagnets. The Ising model is also 
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Beecsds ja model forva lattice fluid of classical particles. 
Recent reviews of the Heisenberg and Ising models are given 
by Rushbrooke et al. (1974) and Domb (1974) respectively. 
The XY model has been extensively reviewed by Betts (1974) 
and much of the description of the XY model in this 

chapter is gleaned from that article. 

In the Sections 2 and 3 following we demonstrate 
the conditions under which the XY model describes an 
insulating magnet and a quantum lattice fluid. Section 4 
describes what it means to dilute a spin exchange system 
with classical impurity atoms and what physical situations 


are described by this Drocess. 


sce Mo Meonetvrc Insulators 

Betts et al. (1970) have shown, using perturbative 
arguments, how in the presence of a crystal field of 
particular symmetry, the isotropic Heisenberg Hamiltonian 
reduces to the XY Hamiltonian. Their argument is as follows. 
Consider a crystal of magnetic ions of high half odd 
integer spin embedded in a lattice of non-magnetic ions. 
The non-magnetic crystalline environment will introduce 
single ion terms in the spin Hamiltonian. For most 
Symmetries (except cubic) the leading single ion anisotropy 
is proportional to (sys . If the magnetic ions are assumed 


to interact according to the isotropic Heisenberg model, then 
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the Hamiltonian for the complete system is, 


SS ACORN AV Se a (872.1) 


neetaiene. Critical region, and for D.o> J > 0 only. the,lowest 
Kramers doublet with S* = + 1/2 will be appreciably 
populated resulting in an effective spin-> System. Ina 


system of two spins for example the XY interaction is 
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O1c-10 
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: ee : 
If the magnetic ion is for instance by? then the ratio 


of the XY non-zero matrix elements to the Ising non-zero 
matrix elements. is: approximately 128:1 . . Clearly then 
magnets which are XY like can conceivably occur in nature. 
In fact some cobalt-halide compounds definitely nave an 
antiferromagnetic transition which is theoretically well 


explained by the XY model. Examples include CoBr,-6H.0 and 
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31 
CoCl,*6H.40 in two dimensions (de Jongh et al. (1974)). 
More recently Algra et ae 1976) have investigated the 
antiferromagnetic transition of the cobalt pyridines 
Co(C.H_-NO)-(C10,), and Co(C.HeNO) §(BFa)5. For these 
compounds the magnetic lattice is simple cubic and the 
specific heat near Us closely fits the series expansion 
results of Betts et al. (1971). 

Theoretical work by Jasnow and Wortis (1968) on 
the anisotropic classical Heisenberg model (obtained by 
faking the limit: as! $>o-ing (3.1.1) )e indicates that if 
a" << gt then tne symmetry of the order parameter is likely 


to be XY like. Universality asserts that this should also 


be the case for the spin-z XY model. 


3.3 Quantum Lattice Fluids 

Yang and Lee (1952) have shown how the spin-5 
Howings MoOdemaconnes ponds: to, a classical. J atti ce, Fuld.» pein 
this context spin up is identified with the occupancy of a 
Wigner-Seitz cell by a particle and spin down represents an 
amoty well... Ihe, lsing Hamiltonian. 1s then essentially the 
potential energy of the system. In a like manner we show 
below how a quantum lattice fluid is represented by the 
anisotropic Heisenberg model (following Matsubara and 


Matsuda (1956)). 
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Fins tiyeeach site. (or cell) is assigned) two 


possible states such that 


ul = ay represents an occupied site 


and V iG) represents a vacant site. 


The creation of a particle at the rth site can be effected 


by the operator at defined by 


- r— in = 
ry Uy and a nUn 0 (Sr oes) 


ee oe T 
The destruction operator a, is the Hermitian adjoint of ans 


and in terms of matrices we can write 


hei AU Ge! re gd 

a, = (6 0? and an (, 9? (3) 23'32)) 
The number of particles in the rth cell is eas aa a 
which has eigenvalues 0 or 1. Creation and annihilation 


operators of different sites commute. To ensure that each 
€e11 1s only singly Occupied, operators on the same site 


obey the anticommutation relations 
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The connection between ay, and ay and the continuum quantum 


field creation and annihilation operators v'(r) and w(r) is 


eee 1/2 
an v (r)v, and a. v(r)v, (3/824) 
and 
& 2 T e 
Drs Wehr) cnavgy Cela) 


where ve is the cell volume. The potential energy between 


two particles at r and r”~ is assumed to be 
Le ial aly AER fe (3.3.5) 


where 6 is the nearest neighbour distance. Hence the 


potential energy operator for the total system is 


The sum in (3.3.6) is over nearest neighbour sites only. 


The quantum field theory kinetic energy operator is 
T= -(A2/2)0 wT (nr) vou) dr (sa2e7) 


where m is the molecular mass. For a discrete lattice the 
integral is replaced by a sum over sites and the Laplacian 
by the finite difference operator and ¥'s by a's. Thus we 


get 
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2 eRe : 
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(fh d/2mé°“ ) ( ibe : (aa,- : aa -q : nn) 
(ron) 
where d is the dimension of the lattice and g the 
coordination number. 
ihesdrandi partition function for the compete 
system is then 
Cru) = Tr exp B{(u-n@d/ms2)N + Ae d/aqms* 3 aula apart a’) 
Gar == OO 
i 
3 z non a} (3 3 9) 
0 none els 
where Ness nh 
r Lf 
‘ xX F ed) ae eee, ny: 
Now if we put ay = ee ~ ish and ay = tp + os ay we can 
make the identification between (3.3.9) and (3.1.1) by 
putting 
am HZ = u + quo/2 - tod/ms 
m= 0» in = ned /qms* 
and ie wae 2 (273010) 
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The pure XY model then corresponds to the neglect of the 


potential energy between particles on adjacent sites. 
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The above derivation is due to Stephenson (1971) 
ana 1s outlined in Betts (1974). Historically: Matsubara 
and Matsuda (1956) were the first to derive the anisotropic 
Heisenberg Hamiltonian as a model for the superfluid-normal 
fluid transition in He’. The XY model of He’ has been 
investigated by Ditzian and Betts (1971) and Betts and 
Lothian (1973) using high temperature series expansions. 
Their results were in good agreement with the available 
experimental data for various critical exponents and 
amplitudes. The results for the critical exponents give 
aco 0. Ocwand, a Vogritnmic divergences Of the= spect fic 
heat (i.e., o~ 0). Dekeyser and Rogiers (1976) confirm 
this estimate of y using high temperature series expansions 


with one more term than those obtained by Betts and 


co-workers. 


3.4 Dilute Models 

By dilution of spin exchange lattice modells such 
Scuthem.\unodel we meadne the intnoductponein the slatticesof 
non-interacting particles. The four common types of 
impurity dilution are classified as quenched or annealed 
systems by their thermal distribution and bond or site 


models by the mode of occupancy of the non-interacting 


particles. 


om at ; 7 \ ¥ Ay ee i r y , a an 
% Be Rah dial, | ; f my ", 7 iv 7 
a ene fe: Hey sion? ge? sob: aH ator eroR: a 7 
Tae . ' re ; sh, ; % ; i) “a ay ; 
a dee PEM ON Fy oer oe # speKeh eles a 


donee d$as: Sar evepah os ae ri she. Ae 


sntoag-Cieliapene. Gey 197 J gbeats ee 


PoAkGea doe 188 Svutsemines meri Bagge 


> 

rf 

— 

2 

= 

it 

@ 

a 
te 
- 
> : aes 
a 


- i oo i 7 sibel = § : : . . a a = 
ST DT Res SM iG, Ta) coed Caper A 1B 
y ad 


rh Pees CRO) avs at Ak awh 
ONE ROSEY A Pe aes ONL IE wSeoe? re ome “yy: 


‘Ba Ava? ott DeeietGo stones nets a4 


ApOR wleben GLITPel Sencniws Abo Teo nol Ee : 
‘ ¥ i 4 7 7 . ae i 
We: aa72051 -o0 1) NEL SCUboge (GRO aoe Ae See 


to <sesnet ndeiea uot et? .pabe RRR) 
pate saig we Matlonpay len, Varth egn!'s| ae mss 
ne tae liorsud tg 2h  bamany Dok: 


7 


| | pibsndinniiny abi rey ie RT PeqaS SO 0. ll - 


y ed , . bade” aa 
ra eee: 


36 


Bond dilution is when the impurity particles are 
placed in the “bond” between two lattice sites. Site 
dilution is achieved by replacing the magnetic ion on a 
given site by a non-interacting particle. In quenched 
Cilutaon the tmpurttves are “frozen in" the lattice ina 
random distribution. On the other hand annealed impurities 
are thermally mobile and assume their equilibrium distribu- 
tion at all temperatures. 

For crystalline magnets it seems unlikely that 
annealing of impurities can be achieved since this would 
imply that the impurity atoms were able to move througna 
the lattice. Quenched impurities in magnetic crystals are 
physically realisable though. The quenched bond case 
could correspond in a super-exchange magnet to the replace- 
ment of the interstitial ion by an ion that does not allow 
exchange. The quenched site case corresponds to the 
replacement of a magnetic ion in a direct exchange magnet 
by a non-interacting particle. 

Annealing of impurities can of course readily 
occur in fluids." Thus it seems reasonable to use the 
annealed site model as a model for a diluted lattice gas. 
There seems to be, however, no mechanism of "exchange" in 
lattice fluids that could correspond to the super-exchange 


mode in magnetic crystals, hence bond dilution seems unlikely 


for Vacti cert fuids. 
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In the next three subsections the site dilution 


problems are considered in detail. 


3.4.1 Percolation theory 

In the next subsection on quenched site dilution 
we make frequent reference to results from percolation 
theory. in this subsection we define the percolation 
problem and its relation to the problem at hand. 

Considenva Vattice with a fraction Dror ats sites 
occupied by magnetic ions. The occupied sites fall into a 
number of clusters and two such sites belong to the same 
cluster if there is a chain of nearest neighbour occupied 
sites connecting them. The mean cluster size S(p) will 
depend on p and one instinctively expects that S(p) will 
increase linearly with p and be a maximum for p = 1. In 
the case of an infinite lattice however, the cluster size 
is no longer bounded. There exists’, in fact, a critical 
percolation concentration P. which is the largest value of 
p for which a given occupied site certainly belongs to a 
enuster of finite Size.. “for p> p. there iS a non-zero 
probability that a given occupied site belongs to an 
unbounded cluster. The relation between the percolation 
problem and the quenched site dilute magnetism is contained 
in the following statement proved by Rushbrooke and Morgan 


* . 
(1961). The limiting concentration p. of magnetic elements, 
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below which there is no critical temperature is the same 
for all such models, whether Heisenberg, XY or Ising and 
regardless of the spin value concerned. This critical 
eoncentration is. of course the critical percolation 
eoncentration.. Win fact “the properties of the 208 
temperature Ising model led Kasteleyn and Fortuin (1969) 
to establish the following analogy between certain percola- 
tion functions and the thermodynamic functions of a 
ferromagnet. 

The pair connectedness (which is the probability 
that two sites belong to the same cluster) is like a 
pare correlationrunction. The percolation probability 
Cthes fraction. of: sites. contained. in. infinite. clusters.) 
has properties similar to the magnetisation. The mean 
Ss tenisize15(p,) is anadtogous to the zero field suscep-= 
Plouiityee Consequently it 1s, possible to assign critica) 


exponents to these percolation functions. For example 
* en ey, 
S(p) = Of(p-pi)/p.d * 


The percolation problem has been extensively 
reviewed by Essam (1972). Table 3.1 gives the critical 


percolation concentrations for some common lattices. 
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3.4.2 Quenched site dilution 

For magnetic Watice models quenched site dilution 
is the occupation of a random distribution of lattice sites 
by non-magnetic impurities. Both the quenched site 
Heisenberg and Ising models have been studied extensively 
by many workers. 

Exact results known to date are few, and apply 
PaiyecOntie Ising. model. The. most importantlof these 
Beoiiceis dle tooGrighiths(1969).) Hissresultusaysechat 
the magnetisation fails to be an analytic function of field 
H at H=0 for a range of temperatures above that at which 
Spontaneous magnetisation first appears. Wortis (1974) 
has observed this occurrence of spurious singularities in 
M for the one dimensional random dilute bond and site Ising 
models. Also Rauh (1975, 1976) has studied a randomly 
dilute spherical model and observed Griffiths singularities. 
Griffiths and Lebowitz (1968) also have derived some exact, 
but véry general” results for the’ Tree energy of the annealed 
site Ising ferromagnet. McCoy and Wu (1968) have studied 
the two dimensional Ising model with a particular distribu- 
tton-of impurities. Behringer™(1957)" was the first=to 
investigate the quenched site Heisenberg (S=5) model and 
found rougly how ie varies with the concentration of magnetic 


ions p. This line of critical points was compared with some 
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eiay experimental results of Fallot (1936, 1937) and 
Forestier (1928) on various alloys with one magnetic 
component. Brout (1959) also worked on this problem and 
Was tne first to point out that Heisenberg like systems can 
be diluted by either annealing or quenching the impurities. 
For recent series expansion studies of the quenched site 
Ising and Heisenberg models, the interested reader is 
referred to recent papers by Rushbrooke (1971), Rapaport 
miav2a.b), Rushbrooke et al. (1972), Eltiot, and Saville 
(1974) and references therein. Because the Heisenberg, i tee 
and Ising models are independent we will not describe in 
detail results of previous work on the Ising or Heisenberg 
quenched site models. Instead, we will quote in later 
chapters only the results needed for comparison with our 


Own. 


3.4.3, Annealed site dilution 

Annealed site impurities are not "frozen-in" as 
are quenched site impurities, but are allowed to assume their 
thermal equilibrium distribution. It is unrealistic to 
expect this et of dilution to be realised in crystalline 
substances. However annealed impurities in a lattice fluid 
could represent one component of a fluid mixture. 

The impurity particles are classical and serve 


only to occupy lattice site which could otherwise be occupied 
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B¥eoneinteracting. particle. .Since the: XY model gives. a 
good description of the Pree aaa de wong Fe UsGest dts, latO 
in He® it is natural to. present the annealed site XY model 


ouiae mixtures. The kinetic mobility and 


as a model for He 
statistics of the He® particles are absent, but for 
temperatures well above the tricritical temperature (where 
phase separation first occurs) these may not be important 
considerations. Formally the tricritical point is where 
the phase separation order parameter, X3> and the superfluid 
order parameter ww, simultaneousiy become critical. 

To demonstrate that the annealed site spin-5 XY 
model has tricritical behavior we establish the equivalence 


of this system with the Takagi model of He?-He@ 


mixtures, 
following Reeve (1976). The argument is similar to that 
used by Wortis (1974) to identify the Blume-Capel model 
(Blume (1966), Capel (1966, 1967a,b)) with the spin-z 
annealed site Ising model. 


Consider the model with nearest neighbour 


Hamiltonian 


(3 7..4551:) 


proposed and solved in the mean field approximation by 


Takagi (1972). Total spin, S=1 for this model. The 
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operators a; and a. raise and lower, respectively, the 
value of the spin projection on the ith site between the 
values S* = +1 and $7 = 0 only. 

Dekeyser and Rogiers (1975) showed that in the 
himit«as Hy - Ho al Ae is equivalent to the spin-5 XY 
model given by (1.1), with mH = Hy 7 Ho . However, since 
the ait and a, act only in a two dimensional Subspace and 


the $3 are diagonal, we can decompose aL into the direct 


Sum of two operators. One of these operators is the spin-> 


XY Hamiltonian uaa given by 
m 
HAY = -d £ (aya,ta,ac) - mH z $2 
m <A J J j 


with mH = Hy + Ho > Where Noa is the number of particles in 
the system. The second operator is the one dimensional 

quantity NCHo-H, )- Since these operators act in different 
subspaces of the spin-1 system they:commute. Hence we may 


write the partition function for the Takagi model as 


N 
Zz = x exp{@N (H-H,)} Tr exp{-BHy/3 fend) 


where ni is "the total number of lattice sites.  Identifica- 
tron of wu. the chemical potential per site, with H-H, shows 
that (3.4.2) is precisely the grand partition function of 


for the spin-z annealed site XY problem. 
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The Takagi model has been extensively studied, 
via high temperature series expansions, by Dekeyser and 
Rogiers (1975) and Rogiers, Dekeyser and Quisthoudt (1975). 
Their results demonstrate, as do mean field calculations, 
the existence of tricritical behavior qualitatively 


3_uet mixtures. 


enaracteristic of He 

Other thermally dilute systems which have been 
studied are the Syozi (1965) model (Syozi and Miyazima (1966)) 
and the bond annealed spin-> Ising model (Rapaport (1972) 


and Cox’ et al. (1976)). 


3.5 Dilution and Critical Phenomena 
In this section we relate the general theories 
given in Chapter [I to the problem of dilute lattice models. 
The effect of impurities imposed in the XY lattice 
model will be to lower the critical temperature rougly 
linearly with increasing impurity concentration. A 
qualitative plot of T.(q) versus q, the impurity concentra- 


Pune cuss diven.in fhigure 3k. 


Seo. eda no! Tneory 


The scaling theory outlined in Section 2.3 
predicts that if we approach an arbitrary point P from the 
high q direction or the high temperature direction the same 


set of exponents should be observed. This is because ordinary 
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Y 


Schematic diagram of temperature T versus 
impurity concentration q for diluted spin- 
exchange models. The arrow A shows the 
direction of high temperature approach to the 
arbitrary point P and the arrow B shows the 
high q direction. 
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scaling theory can be applied to the line of critical 
points, simply by inserting the density (or fugacity) 


Varlabwe as a non scaling parameter in (2.3.1); 


S.o.c. Fisher renormalisation 
Fisher renormalisation of exponents applies to the 
annealed site case and will be practically unobservable 


Since a= 0 


Seoevoy Universality 

According to the universality principle all points 
along the line of critical points are equivalent. This 
means that the critical exponents should not vary with q. 
Weak universality on the other hand demands that ratios 


like y/v be invarient with respect to q. 
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CHAPTER IV 
THES EEWITES CLUSTER EXPANSION 


4.1 Introduction 

Although the first generally applicable expression 
OT ithe finite cluster theorem is due to Domb (1960), the 
method of Rushbrooke and Scoins (1955), in turn derived 
from an application by Fuchs (1942) of cluster integral 
theory to the Ising problem, is substantially similar. 

Nevertheless the power of tne method was not fully 
realised until Domb and Wood (1965) derived the high 
temperature series for the Heisenberg ferromagnet using 
the finite cluster technique. Previous to this work though, 
Rushbrooke and Morgan (1961) and Morgan and Rushbrooke 
(1961, 1963) developed series expansions in inverse 
temperature and density for Ising and Heisenberg ferro- 
magnetis with quenched site impurities. These were 
essentially applications of the finite cluster method, a 
point later elucidated by Rushbrooke (1964). 

The original derivation of the finite cluster 
theorem is reviewed in Domb (1974) and parailels the 
classical cluster integral theory developed by Ursell (1927) 
and Mayer (1939). This derivation is useful in that it 
emphasises the equivalence of the finite cluster theorem 


and the Mayer cluster expansion. However we choose to 
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48 
follow a more elegant proof due to M. F. Sykes and reported 
by Essam (1966). Before we can digesitgetnis “proof gthough a 


few graph theoretic definitions are required. 


4.2 Definitions of Graph Theory 


The definitions presented are outlined in Sykes 
et al. (1966). 


Illustrated below is a graph with seven vertices, 


three edges (bonds) and four components. 


S 
— 


A graph of one component only is cailed connected, otherwise 
it is disconnected. 

Two graphs G and G are isomorphic if there is a 
one-one correspondence between their vertex sets V and V~™ 
such that corresponding vertices are joined by edges in 
one of them only if they are joined in the other. 

A graph Hvis a subgnaph) of Gowheny tnesvertex set 
V(H) is contained in the vertex set V(G) and all the edges 
of H are edges of G. 

A section graph H of 7G hs a subgraph of Gi such 


that its edges are all the edges of G which connect two 


vertices of H. 
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Any subgraph Gp of G which is isomorphic with a 


GGavieG 1S.Said- to, represent a-weak embedding ofegquin G. If 
the graph @ 


is a section graph of G then G represents a 


strong embedding of g on G. The number of subgraphs of G 
isomorphic with g, denoted by (g;G), is called the weak 


Haeerce, constant.of g im G. Similarly the number of; section 


Graphs of G isomorphic with g is’ called the strong lattice 
constant of g on G and is symbolised by [g;G]. 


As an example the subgraph 
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Denoting by, {Ci}, the complete set of connected 
graphs, ordered such that (C53C,) = [C305] + ORT Crait ety. 
results in a complete set of connected graphs in dictionary 


order. The first four elements of such a list must be 


followed by either 


the choice being inconsequential. The final and most 
important definition for our purposes is the following. 
If ®(G) is a function defined for any graph G, 


then 6(G) is extensive if and only if 
O©(G) Ul G7) ,=O(G)e +) O(.G7) ER 


where the graph G U G~ is the graph formed by regarding 


all the components of G and G* as constituting a single 


grapn. 
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In the present context a graph G has. a vertex 
set V(G) which represents the sites of a physical lattice 
having interacting particles ate them, while the edges are 
representative of the interaction between pairs of particles. 
The function $(G) is any thermodynamic function that can be 
defined for the graph. For example it may be the free 
energy ,-Susceptibility,. etc. 

Armed now with these definitions we can proceed 


to demonstrate the finite cluster theorem. 


4.3 The Finite Cluster Theorem 
4.3.1 The theorem 
Let 6(G) be any function defined for a graph G 


which has the extensive property (4.2.1). Then 


where C. € {C.}, I, (G) is either (C. 3G) or [C. 3G] and F(C, ) 


TS a ituUnction specified by the C. Onive 


Aone ive proof 


Since (G) is extensive then from (4.2.1) 


k 
o(G) = z o(C.) 1.(G) (U3) 


where 1 (6) is the number of components of G isomorphic with 
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. a * 
oe and k is any integer such tnat 1, (G) = 4) sfOY wi ake ONY 
connected subgrapiis or section graphs need be considered, 


allowing us to write 


1, (G) = Aw sate CG) 


where Aaj is either (C,5C,) or [Ci 3¢.] depending on the 
enojce Torn the, lattiee constants. . The idictionary sondering 
of {C. } means that Aa; heh Oy dl >a anid Neg = le siov- that 


the determinant of (A...) is one and we can write, 


ii 


" k 
1, (G) = 3 Taped aa(sGe) 


* 
ie Substitution of m5 (G) back into 
(4.3.1) gives the desired result, namely 


k 
o(G) = 2 f(C,) m,(6) (4.3.2) 


with F(C,) = E o(C;) Psi avesis)) 
This finite, cluster theorem holds for beth finite 
and infinite graphs Ge An cinfinite wraph is a ‘physica 
lattice in our case. The lattice constants on a physical 
lattice are denoted by [C;3L] and (C.3L), where L could 
Venuesent ine: f.c.C.,tDec. Ce ance. etc. (Neve dire 


defined as the number of embeddings of C. on L per site of 
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L. Accordingly in the thermodynamic limit the finite 
cluster expansion is 


. kK 
OC = hamit: x 


oo VCC eie GE) (4.3.4) 
Koo j 


As we show in a later section, there are properties 
of the F(C.) which enable us to take k finite and retain an 
exact perturbation expansion for o(L) in certain independent 


variables. 


4.4 Construction of the Expansion 
In practice there are several steps we must follow 

in order to establish the finite cluster expansion for a 

specified function (L). These are tabulated as follows: 

(i) Establish.a dictionary ordered. graph list. 

(ii mb ind the matrix elements Aaj which should usually be 
the strong embeddings, since this enables the list 
in-iu) to be shorter than thagvrequired if thesweak 
embedding scheme is cnosen. 


.) 


(iii) Obtain the weights f(C,;) in terms of the o(C, 


samo (4:33). 
(iv) Calculate the lattice constants M,(L). 


(v) Substitute the F(C.) and 1, (L) it OM (Aino seen) a 
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In order to clarify the discussion of the 
previous section we will derive the expansion up to and 
including clusters of three vertices. The appropriate 


Guapn 1ist is 


and Cy = WAN 


and for convenience we restrict considerations to the strona 


scheme. Thus Ai; = [C;3C,] and m,(L) = [C,;L]. The matrix 
(A)i5 denoted by A(3) is by inspection 
iWeceone 
Oe lie23 
Sei Og OSA 
0001 


which has the inverse 


woe 


it al 


ia) el) 3 

Omeles2 =c 
T(3) = 

Gae04 1) 6 

Doe Or al 


maerequations (4.353) then yields 


f(C,) = o(C)) 


F(C,) = @(C,) - 29(C,) 
F(C,) = o(C,) - 20(C,) + o(C,) 
POC Tec > JOCC ys oCt 7) (4.4.1) 


It is worth noting that we could have arrived at 
(4,451) by. recursivinly solving equation (4.3.2). for 


baa C G = Cos etc. and in practice this method is often 


1 
preferred since it is easily programmed for a computer and 
doesn't require knowledge of the Taj explaicithy. stire 
equations (4.4.1) demonstrate that the F(C;) depend only 


on a(C,) Wiech kK <, ol cand. not.on o(C,) (Oe Me ores Patce . Sieseorabe 
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the essential property which enables perturbation expansions 
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sophisticated (see for instance Martin (1974)). The results 
we need here, for a lattice of co-ordination number gq, are: 
poe = Os) [Ct ean g/2 whe, sk] = q(q=m=1)/ 28 ade PCaeL] = 
nyo.e (nel lattice! Leis anyof the regular Bravais: lattices, 
namely the face centred cubic (f.c.c.)(q=12), body centred 
Expres D Ac. ce )/(q=8) simp lee-cubitcs (5s. c) )'Cqe6)s,. pihane 
triangular (p.t.)(q=6) and the simple quadratic (s.q.)(q=4). 
The parameter m is zero except for the close packed f.c.c. 
(m=4) and the p.t. (m=2) lattices. These data, together 


with equation (4.4.1) are substituted into (4.3.2) to give 


N teh) = o(C,) + ((C,) - 20(C,))a/2 


+ (o(C,) - 20(C,) + 0(C,))q(q-m-1)/2 
s (9(C,) - 36(C,) + 36(C,))qm/6 
gem: OO Mee : C4 ar Zo)e 


4.5 Perturbation Expansions from the Finite Cluster Theorem 
Tne derivation of expansions in magnetic ion 
concentration and inverse temperature for Heisenberg like 
Magnetic models is readily achieved using the finite cluster 
theorem. By way of example we will derive series in quenched 


site density, annealed site fugacity and inverse temperature. 
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Case I: The quenched site problem 
In this problem, applicable to Heisenberg like 

ferromagnets, we add to some of the lattice sites non 
interacting impurities to prevent the occupation of those 
sates by-magnetic ions. The distribution of these impurity 
atoms is assumed to be at random at all temperatures, 
physically representing the quenched defects one might 
conceivably expect in a real crystalline ferromagnet. If 
the probability of magnetic site occupation is p, then 
the system is mathematically represented by weighting the 
factors FCC, ) ane 4 23a) Dy an where Vs is the number of 
vertices of C.. This ensures that all of the vertices of 
Cesare Occupled for f(C.) to contribute in tne expansion, 


‘ 
which 7s 


N-To(L) = po(c,) + po (a(C,) - 20(C,))q/2 


+ p2{(o(C3) - 20(C,) + (C,))(q-m-1)q/2 


+ 


(9(C,) ~ 30(C,) + 30(C,))qm/6} 
ee ee ae (4.5.14 


Usange 4.4.2) directly. 


Clearly (4.5.1) is an exact expansion since C, 


contributes only to terms of order v. and’ ni giie r* ane, 
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Case II: High temperature expansions 

Derivation of series in inverse temperature from 
(4.4.2) is just as simple as for the previous case, but the 
exactness of the expansion is not obvious. For this reason 
we shall demonstrate the method for the specific case of the 
spin-y XY model. 

One quantity naving the desired extensive property 
is the dimensionless free energy Per site f(T) = InQ/NkpT, 
Where Q = Tr exp{-@H} is the canonical partition function 
for the system. For the XY model, the method by which we 
Gotarm 0 for finite clusters “1s contained im Appendix A. 


The results we need here are: 


Q(C,) Secu t ec Casi 
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2 + 2exp(2K) + 4exp(-K) 


iT] 


Q(C,) 


where K = BJ = J/kpT And uv “1s. tnes exchange integra in 
between nearest neighbour pairs of magnetic sites. TInese 


equations give an expansion in K to third order, 
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1nQ(c,) = 1n2 - 


InQ(Cy) = 21n2 + K*/4 + 


1nQ(C,) = 31n2 + Ko/2 +4 


1nQ(C,) = 31n2 + 3K°/4 + K3/a + 


which are the functions o(C.). From these we get via 
SOME TOUS neni E CHEEK Pah FC) Sos ane 
F(C,) = Ko 74, From equation (4.3.2) then we get 


NuvinOv=In2,+ ak°/8 + qmko/24 +... (405.04 


Pauation (4.5.2) is correct only.to second order in K 
because F(C.) produces terms of order ra or higher, 
where he is the number of bonds in C.. To obtain an 
expansion exact to third order would require that all 
graphs of three bonds or less be present in tne dictionary 
grap Viste 

Domb (1965) has written down formulae giving, for 
some different -graph topologies, the term in K at which C. 
LMGeTteICOUCrTDULES 4. hk i is a lower bound for this. 

Obviously we could also expand (4.5.1) in K for the 
spin-> XY model and obtain an expansion in K correct to 


second order. 
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Case III: The annealed site problem 

This, the last case we shall consider, is somewhat 
different from the previous cases and is normally applied to 
lattice fluid problems rather than magnetic crystals. For 
consistencey however we retain the magnetic system 
terminology. 

Just as for the quenched site case we dilute the 
lattice witn non interacting impurities, but allow them to 
assume their equilibrium temperature distribution. We are 
thus required to use the grand canonical ensemble, since 
unlike tne quenched site case the density of magnetic ions 


is itself a thermodynamic function 


‘ea Zo din B/oZ, 


wiere: 2 is the ‘grand-partttion function and the fugacity Zo 
of the magnetic ions is given by exp(u_ 8); with ue the 
chemical potential of the magnetic ions. 

For continuous media the grand partition function 


is; for N particles, given by 


— 


\ 
ey) + ae Wai 0e 
n=0 
wnere O., Z- the Canonical pareltiron tunction Toren panei cles. 


On a lattice of discrete sites finite N corresponds to a 


finite cluster and we may write 
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Z* [9,3¢,] Q(9,) GAS on 


The sum is over all graphs (connected and disconnected) with 
less than or equal to vs vertices. Strong embeddings are 
Useduin (4.5.3) because they do not allow sates which are 
connected on C; to be occupied by sites of 9. which are not , 
a requirement made necessary by our interpretations of the 
graphs. 

The free energy is constructed by forming 
Ine (C, ) fowl Lathe Cc as power series in Lie ASwTONe tie 
quenched site case, if our graph list is complete up to 
endsinmcluding. graphs of M vertices or less), then the series 


] M 


for N In (L) will be exact up to order a 
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(4.5.3) then gives 


= 1 + 2Q(C,) 


seedy 2 Oa ara Kc) 


De eZOCC roa c, )) 9s oneal en ec (oe 
it Sucre az Q(c,) + 2 QUe,) 


have used Q(g,) = Q(C,)Q(C,). 


Proceeding as before and using o(C.) = Ine(C; ) 


gives the weights f(C;) as 


ieee Ole) Zo( OC Wye + 2 (OCC. No he 


F(C3) 3 78 


(Cy) 
up toathi 


nv dne(L) 


pee OC Ga) Tira re (OC. reece Otce 


SEGgeer))el> 2a(chrotes ier OCC, )) 


eRe iy deuSQ( Gr NO GCM) (a Cee I) 


rd order. The resulting expansion is 
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+ Z3{q(q-m-1)/2 ((Q(C,))* = 20(C,)Q(Cy) + Q(C3)) 


+ qm/6(2(Q(C,))*> - 3Q(C,)Q(Cy) + Q(C4)) 


3) ore . Se Be . i ay) a a a ae 
" on : | 


=e ; 4 me + i 7 : 
Re ie Da Sa 
a nye ee *' hee Ny i Bar wet) oe thao Este Bape, 5 
a : te x fs of — + 1. t 
< ; - = 
; w © til, 


™~, 
} 7 | 
b 
7 . ed 
s 
4 has 
z 
y 
‘ * } 
Ci } { 
4 # 
4 ' a ee Sa e 
is 
2 . 
f \ } 
< * 
’ -- 
« it \ ) i 
€ 4 : } 3 
i 4 5 > 
3 i 
S 6 fg é ) ag 
. o . 
7 
< ' a ¥ 
és uy = A ; 
Son ( 
5 
' ve f TY, 2) Taek | ee ol ae 
a i : . d : Ly we | 
er i © 
ae 


63 


Fa GAZI (2( OHS NP 20(c,)Q(c,)) + (Qte, 2/33 


tee ees G'S (ay 524)) 


Although the exactness of the expansion is again 
Poporent, It. iS not a trivial task to verify this: generally. 
Nevertheless the previously mentioned analogy of the finite 
cluster and the classical Mayer expansion indicates that 
a fugacity expansion is not unexpected. The earlier method 
of Rushbrooke and Scions (1955) does give a little clearer 
insight into this property which nas also been demonstrated 
by Essam (1967) for the "hard square" lattice fluid model. 

Series expansions for thermodynamic expectation 
values are extensive and so there is no difficulty in 
obtaining expansions for thermodynamic functions for the 
quenched site and high temperature cases. 


Tentatively let us define 


N 
COs ee #2 nie 0, exp(-BH,)/8, (4.5.5) 


as the grand canonical expectation value for the operator 
0 over N particles, where 0, is the operator for the n 
particles with interaction Hamiltonian Ha: Symbolically we 


can express this as 


<0> 0 exp(-8(H,-un))/Trs, r ; exp (-8(H,-un) ) 
aunt) 
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which is formally equivalent to 


ms] 0 2 ee te 
<O>, = -8 Lim1vt —— {In Tr exp(-8(H_-untA0 
N 10 OX | in,H } p(-B¢ Vee sella 


Now if the operator Sa = es ee ect 0, in the 


n 
exponent 1s additive inethe sense that 


U Sy) = P(S,) " P(S,) 


where Sy and S5 are two independent systems, then the 
quantity In Sea exp(-8P, ) is-extensive since’ the 
commutator [P(S,); P(S,) J = 0. The additive property of P 
holds for all the quantities we consider because 0, can 


UsudieypDe wreitten..in’ tne forms.  0..,.with, operators, not 


ot ae ae | 

Ue) 
physically connected on the lattice not interacting. 
Equation (4:5.5) then is a good choice of 


definition and is the same as 


N abe 
<0) = eae (Ouse (A e526.) 


where 0, is the unnormalised canonical expectation value of 

O over n particles. We note here that the usual definition 
of the grand canonical expectation value <0>, as found in 
standard texts (for instance Huang (1963)) 

is not extensive. This is because the normalised canonical 
expectation value of 0, is usually used in (4.5.6) instead of 


the unnormalised quantity 0, 
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Equation (4.5.5) then is a good choice of 
definition, but is not tne same as is usually used. For a 
graph G the equivalent definition is 

V. ie) 
nes [g,5G]0(g,)/2(G) (4.5.7) 
Sj 
The unnormalised canonical expectation values O(g.) are not 


extensive and add like 
UeGeuy G7) = (0067006) + 0(G7)/0(Gs)) O(G)o(es) 


remembering that Q(G U G*) = Q(G)Q(G"). 
As an example we get for the graph C3, 


320(C,) + BE (TL Gas) Qanl- C3) + 20(C,) + Panues) 
<O00G.)>.= =e 
ve ei OC Ee.) eeenOces i) aez-00e.) 


Expanding <0(C,)> for all the grapns C. in powers of Z 
results in an exact series expansion for N7b<a(L)> as we hac 


before for nv! 


4.6. Inyversionot the Finite. Cluster iheogrem 

In Section 4 it was mentioned that the equation 
(4.3.3) is commonly solved by recursion of (4.3.2). However 
in principle it is possible to invert (4.3,2) analytically fo 
obtain the theorem as 


wv lo(L) = = O(C,)A(C.5L) (4.6.1) 
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This result was established by Rushbrooke (1964) 
for the quenched site problem, but is true in general. One 
other definition is necessary before explaining his result. 

The perimeter s of a graph G- embedded in a graph 
G is the number of vertices of G not themselves in G, but 
connected to at least one vertex of G~. An example is the 


graph C3 embedded into 


as 


which has perimeter s=2 (the sites labelled 1 and 2). The 
number of ways a graph C. can be strongly embedded in a 
graph G, with perimeter s is usually denoted by LC; 3G], oe 


Rushbrookes' result is then 


Vi. 
A(G,3L) = 2 p ‘(1-p)"[¢,sL], s (4.6.2) 
S 

The variable p is the quenched site probability and the 
series is truncated at s or m, where m is the maximum number 
of vertices considered in the initial graph list. The 


result (4.6.2) is valid when p=1, so the finite cluster 
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expansion becomes 


: (4603) 


where S must be greater than or equal to m - eg A Oiaest ie 
term doesn't contribute. 

Counting the perimeter for each lattice embedding 
is more difficult than the normal lattice constant problem. 
However we have chosen to use (4.6.3) since it eliminates 
the tedious recursion which must be performed each time we 
change the temperature in a density or fugacity expansion. 
Tables of strong lattice constants and perimeter counts for 
Bieeee CC. 5D. Cy Ces SrC., DP. tev-and.s.q. lattices. can: be 


found in Appendix B. 


4.7 Advantages of the Method 

The finite cluster theorem is not the only method 
enabling density expansions to be generated, but it is as 
yet the most powerful. 

In the terminology of magnetic systems the 
following advantages can be listed. 
(i) An expansion in density or fugacity contains complete 

temperature information. 

(ii) For the quenched site case an expansion in inverse 


temperature contains complete density information. 
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(i711) When used in conjunction with the XY model all series 
contain complete parallel field information. 
(iv).Using the expansion in the strong scheme necessitates a 


Shorter graph list than most other methods. 


CHAPTER YY 
PP MOUCNCHED USE MODEL 


mei introduction 

In this cnapter we present the results of an 
investigation into the quenched site spin-+ XY model on 
the regular three dimensional Bravdis lattices. 

Lien OD hamwOtestis. Chapter 15 Jas fo ll ows, oln 
Section 2 we discuss first the generation of series 
@expansians=in, p, the concentration of magnetic atoms, 
for the zero field free energy and for the zero field 
fluctuation in the long range order. The coefficients are 
determined for arbitrary temperature. Next we discuss the 
generation of nigh temperature expansions for the same two 
quantities. 

The methods of analysis are outlined in 
section 3. 

In Section 4 we investigate the behavior of 
the critical exponents of the order parameter fluctuation. 
One issue we consider is whetner the exponent, y(p) for 
approach to the second order line at constant p is in fact 
a function of p. We also consider whether the exponent y(T) 


for approach to the second order line at constant T depends 


on T. 


69 


pe 
a 
* 
- 
’ 
’ ~ J 
af 
if 
9, 
a fF 
1 
td 
a ® 
< I. 


' a 
Pals { 
D ; . 
‘ nl 7 a 
: . 1 = 
ar aeul 
a(t a + 
4 , A . js A 
ey - a. 7 Aapirt ‘se 
1 eh 
if i 
' 4 
. : \ 
’ ’ ‘ io < e 
4 ) } 
. ¥ 4 
at 
; i 
iS ” 
‘ 
- ! 
’ 
s 
< si 
*, ‘ 
2 
» 
, “ae * 
i ‘ 
* ¥ / : NO ae 
ap ee 
Ee :  & pion i t et 
z ey 4 eee a ma ii so ae 
Pot aa ‘ 2. ~— 3 y ee ~a' 
_ 0 v7 . ¥ [eG 
i 
~ A di 
‘ + . ke fey a 
Ad * 
‘ Lye \ 4 ees f 


aoa | 
3 
a 
 -_ 
‘ 

re 
2 
— 


& in : 
\ STUL & ; Le? 
: “i 
* " " 
r oth i ¢ i) ye R +" 4 ‘ G 
b’ - * '. 
A 
; { 
e - ——— 
hCa Phe a nant ) ( 
ris as } i ay i 4% 


THanGRS ads seiae fi plies 


>. : . 7 
ahd Tat same, robe Znow 4 OE. i. a 


Py Wee 
: a . ay.) a 7 i a ” i 
7 


L oases hey vain? a: 
: } or vy aly | aa 7 
a ae ne ~ t; = o-. [om se ) 
ad PPR 4 Lt bg od os oiit Te an9 
‘ . is 
‘ 


ih. nit 


: 
= 

3 

7 


¥ oe Vae a eth 


® 


af 


- <p 
lr 
rhe 

a) 


70 
section 5 deals with lattice dependent critical 
properties, particularly the location in the T,p plane of 
tnepsecond’order lames and the amplitude’ of the» critical 
Singularity of the order parameter fluctuation. 


SECTAONLYE9E CONtaImss Summary and discussion. 


easeeneration Of Serves 

In the present investigation we have generated two 
types of series, high temperature series and low concentra- 
tion series for each of two quantities, tne zero field free 
energy and the zero field fluctuation in the long range 
order. 

The method of obtaining these series is explained 


inedenai. ihe Chapter IV. The Hamiltonian for the system: ts 


Say Ste (aja, aca) (seeean 


Aa ei eee Ved 
and the method by which we diagonalise tne Hamiltonian for 
finite clusters is given in the appendix. In Equation 
Cory eK =: BJ 
The expansion for Y¢p,Kk), thes fluctuation in the 
long range ee in terms of p, the concentration of mag- 


netic ions, is of the form 


Yio so © C,(K)p™ (5.2.2) 
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71. 
The c, (Kk) are not of simple analytic form but can be obtained 
numerically to arbitrarily high precision for a given 
temperature. Poa abtes-5.1,5..2 and 523 Sareatabulatedsfor 
a selected set of reduced temperatures, t = K/K, Co through 
Cy for the three cubic lattices. 

The high temperature expansions of the free energy 

moemene t.c.c, (fF), obec... CB) tands:..cy (US) Mattiices are 


respectively 


f-(p 5K) = Wa a" lz Ke + 2p°K> + (-15 p? dh 13 o° 1 4a oy K4 
+ (=p? - 5p) + 105 p?)K? + (sag po + as p> 
- 232 pf _ 17] p? + 308, po)K® + (qh p? + 3aq pe 
- gtd, p> - 6ip® + 96% p’)k’ + 
Eb ia2e a) 
fa(p.k) = In2 + poK’ + (-35 pe - is p> + 5 Bal Ke 


AGG ban DOW 9 1 beOe 2250 P Bae 
(512.4) 
42 Lis eile wh le ee eed 
f.(psk) =In2 + pK” +A-ay PL - ge HBP IK 
To: 85 Le Sah aI ed eS SETS con 
aoe ah {cee See 


omens) 


i ee 
= ; 
; I 4 i. ee yak | RLY, on 
Sa ES Ss: Ag ge Cte 
es te “ 
i‘ ney ia! Break 
4 a ae 4. 5] 
ay be c hubs, nk 


ha ‘ * a 
) 1g 
' > Ta ba : ‘e4 
; ‘a 
= ae ob at ae 
‘ ie Bl 
¥ 
” 
ah ae 4 San 
r y 4. amy 
~4 + e A “= * 
hn, 4 
Me 
‘ i 
F uy 
‘a 
¥ ‘ 
4 4y ’ 
ne es fe th 
; Se | 
e 
Wh; SN 
an + f a 
y rt “4 oy j 
; Vi ws 


Caos ce 


| 


Ue 


9°Z2vEGe- OL OEEs 020° 19S- 

v6OLEL 688°S19 6L0°8El 
292 °0St 250°S8l G620°LL 
Gre "HSI 65195" €L 6EGE HE 
9€€h 6S GEOL *ZE 2618" LL 
Z68E "SZ 2vl6°SI 9176 °6 
$900 ‘ZL £005°8 9296'S 
9/00°9 G0SL*b LSZZ"€ 
06L2°E £8982 pely'2 
SEs" L EVIL L L9/9°L 
ie 9, = 


OLeee0 
pouzuad a9e4¥ BY UOJ (Z°2°G) Ul (¥)% SPUSLILJJIO) 


OE CnLIS 
StS pe) 
€860°62 
0686 SL 


Beets) 15) 
SvOle 9 
V6¢0t P 
v988°¢ 
SAS Nard 
AS b 


(1) 


*99192e1 I1qno 


c6c¢ ¢ SLCOR L°O 
BOG ot L8¢0°9 ¢ 0 
0896 °6 Leee Pp Ean6 
DICE Oma OSE Gc v°0 
06¢26°L £609 ¢ ooo 
CLE VEST? oo 
GO6L°¢ poss l ZO 
v6sl°¢ OSV ISL 8°0 
GEL 93977) 6°0 
OvSr'l 90 ccs 0°L 

ae} 25, 1 

17g arqey 


13 


e622 Wee 156°1S9 vLS0' bP €2S°19 0601°8 0919°S Zo 

e6LbZl 6€8h°S9 GLSE°82 Bilis hese 9/b2°€ 1°0 

Ov8z°E2 9/02°S1 £981 "6 00€0°9 Lose*¢ Loez‘2 9°0 

9518°S 1569'b Z009°€ 2406 °2 eleccd 1969" 1 8'0 

fecaot Geis | 9959°1 €009°L SyLyl 6r9e'l OL 

hs 9. 5 v5 25 ae ; 
pre'O = (L)>¥ Buryea 29193e1 9Lqno 


q 
pauzuas Apoq ayy sos (Z2°Z°G) UL SqUaLoIyJa09 ZG LGR] 


74 


SGreee (L)0 Buryez ao13qreL aLqno 
atdwis ay vos (Z2°2'°G) UL SqUaLdIyya0) €°G SLqR] 


¢°0 


Pa) 


90 


8°0 


vL°620L- v9" 2S GpyL*62- 66€L 6h E8bG°2 7890°S 
v68L°2L 2006 ° 64 se/8'sl E90E "EL 9282"p BLOE*E 
9/00°6L GeOUeel ee toot L L68¢°¢ 26L6°2 SEve 2 
ZOLb'S v9Eb'D pren'€ GLvs'2 VLLO*2 €66/°1 
Oves'l 6y8Z "1 9799" 1 icon | 2GLb"L vSSr'l 
Ly 95 Sy v5 E5 a 


7. 


¥ 
a ees eo ee ine gel 


rr x 


F ena f ; 


ia 7 ey. } ‘ 
Mw ee 
MS kh A te dels: 
Wy, oe 


- 
v ve 
‘ sy ‘ ¢ 
A 
Fg » 
Pw eo 
rs 
n 
- 
a 
i 


4 - 
cv) ‘ 
1... io Ge 
i« 4 


ma 


= = s 


——- = wl 


ao eit 


* 


parameter 


¥e(psk)/p 


Ye (psK)/p 


tes 


The high temperature expansions of the order 


fluctuation on the same three lattices are 
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5.3 Methods of Series Analysis 

Given the first n terms of the power series 
representation of a function F(Z) we wish to determine the 
radius of convergence Ze of that power series and the 
asymptotic form of F(Z) as 4 > i . The two most common 
methods employed in critical phenomena are the ratio method 


and the method of Padé approximants. 


Selene ratioumetnod 


Suppose that the coefficients in the series 


anes (5.301 


are known to degree n and are of the same sign. Assuming 
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tiat all. the a, are of the same sign, the sequence of ratios 
ia apy 1/ay must approach the radius of convergence ae of 
F(Z) as 2+» . If we further assume that F(Z) has-the 


asymptotic form 


E(w el) - 2) aie (ogee. 


as Z>L. then the ratios must have the asymptotic form 
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hiceradius of, convergence of (5.3.1) is the first: singularity 
from the origin in tne complex Z plane. A plot of r, versus 


as thes intercept. onmetne 1/2) = 0 axis:and 


1/29 will yield Ta 
the slope will give z7'(g-1), provided thes form Be5ms.'c)las 
eorrect. By expanding (5.3.2) as a series inv Z the asymptotic 


form of the ay for large 2 is determined as 


a, v Ag(gtl) ..- (gt2-1)/2!Z0 (5.3.4) 


from which the amplitude A can be obtained using the previous 


estimates for g and Z : 


5.3.2 Padé approximants 

The ratio method is usually used to give preliminary 
estimates for the position of the critical point and the 
exponent g. A more widely used method, especially for series 


of moderate length, is the metnod of Padé approximants. 
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A [N/D] Padé approximant is the ratio of two 


polynomials of degree N and D, so that 


| ny . nyZ + es tte 1h 
[LN/D] = GT a ee (5.3.5) 
neds? sod ia Ree badew 


The N+D+1 unknown coefficients of equation (5.3.5) are 
ecuated=to the first Nt0tl coefficients of the series (5.3.1). 
If the first L terms of F(Z) are known, then we must have 
Nt+D+]1 <¢ L. Since the singularities of the Padé approximant 
consists of simple poles only we should represent by Padé 
approximants functions which have simple poles only. If 

the function to be analysed has the simple form (5.3.2) then 
the functions [F(Z)]!/9 and dinF(Z)/dzZ = [F(Z)]~'aF(Z)/az 
have simple poles. If g is known then poles of Padé 
approximants to Maga © give estimates of 1s and the 
residues of the poles closest to Ze give estimates of Awe, 
If g is not known then the location ‘of poles of Padé 
approximants to dinF(Z)/dZ give estimates of Ay and the 


corresponding residues give estimates of g. 


5.4 Critical Exponent Behavior 

At the phase boundary the fluctuation in the long 
range order (or equivalently the initial susceptibility) of 
the spin-t XY model is expected to diverge strongly. For the 


2 
purposes of series analysis we assume a Singularity of the 


Simple form 
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¥(psk) % C(p)E1 - K/Kp(p)]7¥ (PI (5.4.1) 


for approach along a path of constant p in the K-p plane. 
For approach along a line of constant temperature the 


Singularity is assumed to have the form 


V(p ke OCKYET = p/perk) fy (5.4.2) 


Admittedly more complicated behavior cannot be ruled out. 
Indeed various more complicated types of phase transition 
behavior have been established (McCoy and Wu (1968), 
Griffiths (1969)) or made plausible (Suzuki (1974), Harris 
(1974)) for other models. 

Most previous analysis has concentrated on the high 
temperature series, for which the behavior (5.4.1) is assumed 
(Rushbrooke et al. (1972), Rapaport (1972a) and references 
therein). Accordingly, we consider the high temperature 
series first. As critical exponents are lattice independent 
and as high temperature series generally behave best on the 
foc. Chelattice we concentrate on thatlattice in attempting 
to determine y(p) and y(K). QOur-results on the’b.c.c.s and 
s.c. lattice are consistent with but less precise than those 
Onwthertac.c. lattice. 

Detailed analysis of a 9 term high temperature 


series for the fluctuation of the pure XY model on the f.c.c. 
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Batuaece yieldedyy(1),= f733s+ 0402, and Ko (1) FeO uce LO +e OR 000 
Moetts ce. al. (1970), Betts y11974)). Naturally the analysis 
of the seven term high temperature series for Y(p,k) for the 
dilute XY magnet will yield less precise results. 

Assuming the validity of (5.4.1) estimates of y(p) 
can be obtained from Padé approximants to the logarithmic 
derivative of the high temperature series for Y(p,K). The 
results of such analysis for a selection of values of p for 
the high degree, central approximants are displayed in 
Table 5.4, 

The most noticeable feature of Table 5.4 is the 
increasing scatter of the y(p) estimates as p decreases. 

PhS scatter, even -at.p = 1, is sufficientlyegreat to 

indicate that the hign temperature series have not yet 
Eetttededown to their asymptotic behaviors Tf, nevertneless, 
some. “best" estimate of y(p) be extracted from Table 5.4 

and plotted against p then the resulting curve falls roughly 
halfway between the similar curves of Rushbrooke et al. (1972) 
for the spin-5 Ising and Heisenberg models. Given that y(p) 
increase as p decreases for all these models, then the 
behavior XY model is as expected from simple interpolation. 
However we do not believe that the Padé approximant evidence, 


for the XY model, is sufficient to establish firmly a 


dependence of y on Pp. 
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Table 5.4° Estimates of the critical exponent; y(p),° from 
LN/D] Padé approximants to series for 


di Vogeyterki7dkton thetticrc: Tattice 
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Ratio plots also indicate that the high temperature 
series for Y(p,K) have not yet reached their asymptotic form. 

The low concentration expansions of Y{p,K) tell a 
rather different story. We have been able to analyse these 
series by both ratio and Padé approximant methods for a set 
Of Values of t = T/T? (1). Fore 0.7. < tos )S0uthe sata or ol ous 
q@ré.remarkably linear. figure 5.1 illustrates the behavior 
Ometne satios for -a typically good series, at t = 0.9 
Table 5.5 summarizes the results of the ratio analysis. The 
assigned errors are confidence limits obtained by examining 
the range of plausible straight line fits to the points on 
tne ratio plot. 

Alternative estimates of y(t) are obtained from the 
residues of Padé approximants to low concentration series 
eepansions of Y(p;K).on the f.c.c. Vattice. The nesults of 
such analysis for degree 4 and 5 approximants is displayed 
in Table 5.6. Confidence limits’ could be assigned on th'e 
basis of the degree of consistency among different Padé 
approximants to the same series. Note that the central 
estimates drift slowly upward as t decreases. Ihe estimates 
from high denominator degree approximants drop while those 
from high numerator degree rise snarply as t decreases. 

The ratio estimates of y(t) and the estimates from 


central Padé approximants agree well with one another. 


ee a 
er a ; ha on rotate: ae ee *, ; oh bf ; ae 
aoe quia. a O ane oa ‘ahs ae see 
ra dina Lt aby bem es ae ie 


1 [ ath o } G ; i a re ce a aa. ae Fs te ai 

Coeaiuh Soy bs eee eew mete + ned wey 

as ails tah | id q 
j 4 went ig Ha | ‘aunt Seb Aas a iM ve 

ae Sty | Suh. a meh Hoge ga’ ee 

5 | Wings shy a, a 

Av Ro paieta i ; ya 3 me eet? 


=F) 
1 
—— 


Fog at as cae me! . 

dd apaa eines pr ven ee 
caivle netaERe adage wes ‘#2, semi 
oF fu cay SHE fepitisl 30, Livia 
welagtt ah nei Kus 176 id Nagy. or 
ails? 1e boRei eet Sh ig ds: 1 svi BS 


Bick newer Brig une ) 
(sytheaeie Beas S20} eieg ing 
gen ites ant Papdescd ens 2 pinged i 
gents sted an Ae gem xongae aay 

, ae Hen ieb & fe atimiede aid fay :t 

mot; 242 ett ee Fis ‘bre iets dhe 


ASHSONS Bio. ne ku Wie vente ed hoi 


83 


1,20 


1,14 


re 12 


1,10 , 
1/3 1/4 1/Sip te 0 


/n 


Figure 5.1 Ratio of coefficients c /c._, in the low p 


expansion of the fluctuation Y(p,k) versus 
/ne forthe 0422.0) Ks. 059: on ethers ceice 
lattice. 
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Table 5.6 Estimates of the critical exponent, y(t), from 
[N/D] Padé approximants to series for 
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Naively one would conclude that for the undiluted spin-y XY 
model y = 1.23 + 0.01 while as the critical temperature 
drops due*to dilution the critical exponent y(t) increases 
eomrececiwmanvalue: ys) 04 #80) Peat*t™=s0: 5° {orvathoe = 06). 

Both features are somewhat unsettling.: For the 
Biipen ow = lf ce XY model sy e="1. 33°%+40..02" (and= for ther o> << 
emmnoee)| y=" 1.31 +°0 02, Ferer and Wortis’(1972))" “The 
equality of y and y for the undiluted model seems rather 
firmly excluded. According to prevailing views of critical 
behavior, near a line of second order transitions (Griffiths 
and Wheeler (1970), Griffiths (1973)) one would expect the 
free energy and thus Y(p,k) to be a generalized homogeneous 
function of an appropriate pair of independent variables, 
bak) and: g(p,k)s A special direction. in the pk plane is 
@efined by the tangent to the second order transition iine, 
K = K.(p). The X} = 0 axis then must be chosen tangential to 
the second order line at the point of interest (here, at 
p=1) and the g = 0 axis is then in any direction making a 
finite angle with the tangent. If Y(X’,g) were a homogeneous 
function of its arguments then y(1) and y(1) would have to be 
equal, contrary to our findings. 

Secondly, ideas of universality or smoothness 
(Griffiths (1970), Kadanoff (1971), Bettswet al. (197i), 


Ferer and Wortis (1972)) would require y(t) to be independent 


87 
of p. The estimates of y(t), particularly from the ratio 
plots, cannot readily be reconciled with a constant jy. 
Smoothness would also require y(p) to be independent of p, 
but our high temperature series results are not precise 
enough to rule out a constant y 

The above results are all-based on the f.c.c. 
lattice. However the corresponding results for the b.c.c. 
anaeuseC. lattices, though less precise, are in accord with 


Pres rc.c. results. 


3.0 © Lattice Dependent Critical Properties 

From the series we have generated, if they be long 
enough, we can derive a number of critical properties which 
are lattice dependent. In principle we can estimate for the 
Prcec.. Ducac. rand: siic.: Wattiices the famp li tudes: (C @p)vand 
C(t) as defined by (5.4.1) and (5.4.2) respectively, 
analogous amplitudes A(p) and A(t) for the specific heat, 
and the location of the second order line, p(t) or t (Pp). 

To locate the second iorder iine’, p.(t), we use 
the more regular low concentration expansions for Y(p,K) 
rather than the high temperature expansions. If y(t) were 
given the most precise series estimates of p(t) would be 
obtained from Padé approximants to yl /Y, Now in Section 4 
we have noted some apparent t dependence of tne estimates of 


y . However a 10% error in y yields only approximately a 1% 
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error in Pos which allows us to determine p(t) with 
reasonable precision. We adopt a uniform value of y(t)=4/3. 
This represents a reasonable mean value of y extracted 

from Tables 5.4 and 5.5, and also this is very close to 

the value of y(1) estimated from high temperature series. 

Over smosteofethe range 0.3 f<4t".< JF. 0%estimates 
of p(t) from higher degree central Padé approximants to 
y3/4 are well converged, and we have taken as the best 
estimate of Po the averagemo fh giiey! 1/5 lee | 2Y 450 e [S/ 3a. wd ea. 
mee 2/3), 03/2] and? (2.2 approximants. Our confidence 
mimits. given y-4/3," arevagfraction of a percent over the 
above range. The quality of the series deteriorates rapidly 
for’ t < 0.3 however. 

Our results for p(t) are displayed in Figure 5.2 
for all three lattices. For comparison we have also plotted 
the mean field result, pare t, the:results for the S=1/2 
ising mode)’ on the f.c.c. lattice (Rapaport (1972a)) and the 
results for the S=1/2 Heisenberg model on the f.c.c. lattice 
(Rushbrooke et al. (1972)). Rapaport was able to extend the 
Ising model curve down to t=0 while Rushbrooke et al. were 
unable to extend the Heisenberg model results below t=0.4 
For the f.c.c. and b.c.c. lattices we have been able to 
estimate PLtiforstta 0.2, and as the curves are so smooth we 


feel confident in extrapolating them to t=0. Also plotted in 
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Critical temperature, t. = TCp)/T (0) versus 
p for the mean field approximation (M), the 


spin-z Ising* model on’ thee fic. cis lattice: (1)% 


the spin-y Heisenberg model on the f.c.c. 
Tattice.(H) ;.and.the spin-7 XY model on the 


fcc .6n CE), DoCeGa( Bia and: Sec. VS) wacticese 
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Figure 5.2 are the critical probabilities for the site 
percolation problem (Essam (1972) and references therein). 

| The p(t) curve for the Ising model, as already 
noted by Rapaport, appears to intersect the DiaxX Seat tne 
Site percolation critical value, pag with infinite slope. 
In contrast the XY model curves appear to intersect the p 
axis with finite slope. Of course for any lattice model 
of dilute magnetism we must have p.(0) > Pe because there 
Brean. inninite clusters for p.< pe: For the Ising and 
the XY models the inequality seems to be an equality. The 
finiteness of dt/dp at t=0 for the XY model as opposed to 
the infinite slope for the Ising model reflects the 
difference in character of the elementary excitations. It 
requires a finite energy to overturn spins and thence reduce 
the Ising model magnetisation while the XY model presumably 
has spin wave like excitations of vanishing energy. - 

For the Heisenberg model curve to intersect the p 
axis at Po (Ff ttdo SOjmaitamust Nave ja, range of wpward 
curvature, as discussed by Morgan and Rushbrooke (1963). 

In contrast the Ising and XY models seem to have downward 
curvature only throughout their whole range. 

Finally we present in Table 5.7 our best estimates 
of the terminal gradients at both ends of the p(t) curves 
for the XY models and for comparison the figures for the S=1/2 


Ising and Heisenberg models (Rushbrooke et al. CTO7 Zee 
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Residues of Padé approximants to y3/4 yield 
critical amplitudes C(K) or C(t), assuming the validity of 
(3.2) and continuing to adopt the universal value vy CE) = 4/13. 
The results are of course not so well converged as those for 
p(t) from the poles of the same Padé approximants. Given 
the foregoing assumptions we obtain estimates of the 
amplitudes Cy(t) (X°= boB cor S) withea typical econtidence 
iembt OCi- 2 OY 3% over the range 0.5 < t «< 1.0 

Figure 5.3 exhibjtssC(t) for the three lattices. | A 
somewhat unusual feature is that all three curves intersect 
at t = 0.8 and C = 0.95 . We do not however attach any 
particular significance to this feature. 

Attempts to analyse the specific heat series have 
not been successful. This is not surprising in view of the 
shortness of the series and the well known intractability 


of specific heat series in general. 


5.6 Discussion and Summary 

In Section 2 we have discussed the generation of 
sepieseexpansions. for; f, the zero field free enengy wand, Y, 
the zero field fluctuation in the long range order or 
transverse magnetisation of the spin-y XY model with quenched 
site dilution. Low magnetic ion concentration series of 
degree seven in concentration, p, have been generated for both 


f(p,K) and Y(p.k) for arbitrary K, the inverse temperature 
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Figure 5.3 


The amplitude C(t) of the order parameter 
fluctuation’ versus’ t =31/1- (1) sor sche 
fcc. CF)... C..C. Cb) eamide Sacra Se aticase 
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variable and coefficients for Y for selected values of K 

are tabulated in Section 2. Using the low concentration 
series high temperature series of degree seven in K have 
meemegenerated for f. and Y for arbitrary p2. The resulting 
sewmlesitarewliisted in Section ,'2 also. sain Sectionie4 «above 
we have reported the results of analysis of the high 
temperature series for Y(p,K), which are not conclusive but 
point toward a p dependent y . The low p series for Y(p,K) 
are however more regular and yield an apparent p dependence 
of y as illustrated in figure 5.2. 

Assuming a "best" value of y = 1.33 we have 
determined the second order phase transition line 9 (t) from 
tie lowep series: for: Y¥(pykK) for all three. lattices. » The 
results are displayed in figure 5.3 with for comparison the 
p(t) curves for the Ising and Heisenberg models. One 
interesting feature of the XY model curves, in common with 
the Ising model curve, is the downward curvature everywhere. 
Another interesting feature is that the XY curves seem to 
intersect the p axis at the appropriate percolation proba- 
bilities. In contrast with the Ising model, the XY curves 
fhavestinite slope at. 1 = 0. It 1s not possiblesto sayang “he 
Heisenberg model curve intersects the p axis at the percolation 


probability, but if it does then it must have a region of 


upward curvature. 
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The linear behavior of TO (p) as a cUNnCtIONGOF ip 
near tne percolation concentration P. nas been verified by 
Stauffer (1976). He has used the Spin-wave argument of 
Shender and Shklovskii (1975), together with de Gennes' (1976) 
scaling law for the conductivity of a random resistor network, 
to arrive at the relation TA (p) a (p - p.) as p > P : 

Scaling theory for dilute magnets in the percolation 
limit has recently been investigated by Essam et al. (1976) 
ona stautfer: (1975)<9. In, particular: the: results .of Stauffer 
obtained using cross-over considerations similar to those 
used in Section 2.3, predict the divergence of the amplitudes 
mit) as t > 0-shown in figure 5.3.. Finally we, note. that 
Cox et al. (1976) have verified that. for the quenched bond 
Ising model the ratio y/v doesn't depend on p, the concentra- 
tion of allowed bonds, whereas y does. This would seem to 


verify the weak universality hypothesis of Suzuki 1 Naabhsls 


case. 
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CHAPEER V1 
THE ANNEALED SITE PROBLEM 


oO. introduction 
The nearest neighbour spin-5 X¥ model in a parallel 


magnetic field H, is defined by the interaction Hamiltonian, 


H = -J 2 (ara, af a.at) - mH & Sa 
<1) > J J j 


C623) 


The system we consider in this chapter is governed by the 
Hamiltonian (6.1.1), but the lattice contains impurity 
atoms which are thermally mobile. The possible existence 
Seeds cricritical point. for “this model Ws discussed in 
Seton: 3,-3.4 

ines plan of this chapter as as, follows. ginwsection 
2 we tabulate some of the series expansions used in the | 
pee ar Section 3 contains the analysis of the series 
for the three cubic Bravais lattices when H=0 and for the 
fcc. «lattice for some non-zero yalues oft. In Sect ionn4 
we relate the spin-z annealed site XY model to He?-He® 
mixtures and show how the A-lines may be obtained in terms 
of experimental variables. Attempts to locate the tri- 


Gratical point for the f.c.c. .tattice when N=Osanesoutiinedain 


Section 5. A short discussion and summary are contained 


Tm Section 6. 
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6.2 Generation of Series 

The series for the fluctuation in order parameter 
and for the grand partition function are generated by the 
method outlined in Section 4.5. Series in magnetic site 
fugacity Z may be converted to series in magnetic ion 


density va by using the formula 
N_(K,h,z) = ZdIn =(K,h,z)/9z Onc ns) 


where K = 8d and h = mH/J 
Thevexpansion for. ¥(K,h.2).athe #£luctwation in 


long range order, in powers of z when h=0 is of the form 


G6.2.12) 


The coefficients b, (K) are not analytic, but are obtainable 
to anbitrarily «high precision for, a given value of K. 
Tabulated in Table 6.1 are the coefficients Dos PEW indn Men 88) 
tieet.cec. llattice| for selected values of normalised 


temperature t = K (1) /K. 


The series (6.2.2) may be converted toi series 
in magnetic ion concentration n, as discussed above. In 


zero field this may be written as 


XR 
Y = ihe + Ve a, (K) Mp G6.22:8)) 
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The coefficients aos Usennwo/ Tor the» f.c-c,. Vartiucemare 


tabulated in Table 6.2 for selected values of t. 


6.3 Second Order Transition Lines 
In order to determine the locus of critical 


densities ne 


as a function of temperature we rely solely 
on the analysis of the series for the fluctuation in order 
parameter. The only other available function expected to 
diverge at the second order phase boundary (and then only 
weakly) is the specific heat, for which the series we 
obtained failed to yield consistent results. As is usual, 


we assume a divergence for Y(K,h sn, ) near the critical 


density with the form 


Dimon on Vie COR h)ET Spm cek snd Sel) ee. 3unn) 


where the path of approach in the Kn. plane is along lines 
of constant K and h, with h = mH/J. 

Both the ratio method and Padé approximant 
techniques have been used to locate ne for fixed values of 
t = pele = 1)/K. In general the results of the analysis 
are very similar to those presented in I for the quenched 
Sitescase. ‘In particular, ratio-.analysis and the residues 
of Padé approximants to din¥(K,h,n_)/dn_, Poretheretrocrice 


lattice show that the apparent value of y(K,h=0) varies from 
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ere ge 1.21 + .01) fort ='1.0 to y = 1.5 #001 etorets = 0.5. 
Values of y from some fourth and fifth order Padé approximants 
to din¥(K,h,n,)/dn, are presented in Table 6.3. Estimates 
Gravrare scattered for ™% <*0.4°%: 

The value of y({t=1) agrees within the error range 
with the value 1.23 + .01 found for the quenched site case. 
Qn the other hand, high temperature series analysis for the 
puuemcase Gave y = 1.333. h).002 (Betts, Elliott: andeleen( 1970). 
Dekeyser and Rogiers (1975)). According to the same arguments 
from scaling theory presented in Chapter 5, this value of y 
is expected to equal y(t=1) for the annealed site case. If 
Mindeed 7(t=1) is less than y, as our estimate indicates, the 
reason is not apparent. Fisher (1968) has also argued that 
the observed value of y(t=1) will be greater than y if a is 
non-zero zero, regardless of the sign of a. Also if the 
specific heat has a logrithmic divergence at tne second order 
phase boundary then tne observed value of y(t=1) should 
equa ly. 

The variation of y with t was also apparent for the 
quenched site problem. The ramifications of this result 
with regard to universality were discussed in Chapter 6. 
For the quenched site case nowever, the magnetic ion 
concentration is not a thermodynamic variable, so any t 


dependence of y is not necessarily a violation of universality. 
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Table 6.3 Estimates of the critical exponent y(K) from [N/D] 
Padé approximants to series for din¥(K,h,n,)/dn_ 


Gn the focrcasliattice for i=0..0. t=K (ns) /K 


However, for the annealed site problem the magnetic ion 
concentration is a thermodynamic variable and any variation 
of y with t does imply a violation of universality. 

Pore tire, D<CiGerand S<c. Jattuces im zerontuerd and 
porsine f.c.c./ lattice*in non-zero field.,“the estimates of 
y(K,h) are, as expected, not so consistent as for the f.c.c. 
Maccriee in zero field. Table 6.4 contains the values. of ¥(h) 
from Padé approximants to dinY(K,h,n)/dn_ TOY tte On Guan 
Panwous Values of h for the f.c.c. lattice. ‘The® various 
Padé approximants snow no apparent variation of y with h for 
h < 0.4 . For tne largest value of h considered, h = 1.0, 
tne estimates of y are more scattered, but there is stil] 
no significant average deviation from the estimates of y 
at h = 0. Constancy of the high temperature exponent y with 
Varaatton in h for the pure spin-z XY model was also observed 
by Dekeyser and Rogiers (1975). In other words y(K,h) appears 
Bo ebe universal with respect. to H, but not with respect’ to T. 

To locate the second order transition surface 
Bebe), we use estimates of ne from Padé approximants to 
Riciciene yin’, assuming a constant value of y = 4/3 (see 
Chapter 5) for the three dimensional lattices. Estimates 
of ne from the fifth and sixtn order Padé approximants are 
tabulated in Table 6.5 for the f.c.c. lattice with h=0 and 


: @ 
various values of t. In general the estimates of ae from 
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Mapes o.s Estimates of the critical exponent 7(h) from LNy/D] 


Padé approximants to series for dinY¥(kK,h,n,)/dn, 


On, the f.C.Cuslattice for Se eh mAs 
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this method are rather insensitive to small variations in 
the choice of y . This supports our adoption of a universal 
value for y for the purpose of obtaining a transition surface. 
Zero field second order transition lineswin the tne 
plane for the three cubic Bravais lattices are exhibited in 
Figure 6.1. A comparison of these lines with the equivalent 
set for the quenched site case (Figure 5.2) indicates that 
the line of critical density of magnetic sites is of 
opposite curvature to the line for the annealed site 
problem on the same lattice. The terminai gradients of 
these transition curves at t=1 match those for the quenched 


Spce case for each lattice. 


3 “ 


6.4 He~-He’ Mixtures 


Up to now we have interpreted the annealed site 
problem in terms of magnetic variables. However tne problem 
is physically better suited to a lattice quantum fluid, 


4 3 


for example He®-He mixtures in Which= the he; particies) are 


approximated by classical non-interacting particles sesit we 
te 3 4 
let N3 and Ny be the respective densities of He and He 


particles per lattice site and make the identifications 


= zs Ae 
Nn. 1 Ne. (6 1) 


Pee ee See (6.4.2) 


ae 


a bn oi 2 SAE ‘ ta ba | : 
Later ey Eat am nt i eg neeets Muhsin. 
eu a | rs net é ‘oe siialee Neneh, ena 

Ae ety: i ‘ea i Sty Fa aga stele a pe hh : a 
eon dtige ame, “togidoal wiewen® td.) a 

a Beit Ne ie gin be wicoity to aged a ; ay tae 
eiul byt Le , “ewgtt) aay sate) oe sth 4 
: aagre ai Lamune a anyon edt thnk 7 
. nef mie shee seg ant si} 93 wae at : e 
RUB UNE! | ahaa at ‘: -gathtet see ai 


_ 
nab are DN ‘SeoHe ean Hey +6 so0ws aotst ie " 


re on owg 


Ale 
x, 


Poars (Brig vaugiah oni a § uf pidahdiild 
hy FUT: oe foateiet-t be: 19 
Bibs oe me vot rapa & 


107 


0:8 


O26 


Figure 6.1 Critical temperature t = Retna) versus n_ 
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we can obtain the density per site n = nN + Ny and the He® 
concentration X3 = n/n in terms of our natural variable Nia 
The choice of sign in (5.2) arises from the observation that 
<S_(H)> = - <S(-H)>, while the Hamiltonian (6.1.1) remains 
invarient under change of sign of magnetic field. Provided 
that the expansion of <S,(K,h)> as a series in a is well 
converged for each value of K, h and ni, on these ridiicea acviine 
we can find the A-lines as contours of constant n in the tx. 
plane. In practice we have been able to do tnis only for the 
Pe c.ce WattLice using data from the ¥(Khyn) series for 
Becta. enabling us: to find part of the A=Linewwor 
Wee + .02 . In figure 6.2 thessecond-ordermtyansi t10n 
lines in the tx, piane- for h=0, +1.0 and -1.0 are displayed 
TOY the f.c.c. lattice, while the dashed curve shows» the 
Peiinertor n= 0.65 . The lines of constant. h in Figure 642 
have confidence limits in the X3 variable of 2-3%, while the 
constant n curve has confidence limits of about 6% in each 
ordinate value. 

The series we have obtained for Y(K,hyn) in powers 
of n, are generally well behaved for 0.4 < t < 1.0 and 
Mibu 10 » Onithe othen hand the series fons CK np2)esin 


powers of z can be analysed consistently in the region 
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Critical temperature t = K.(x3=0)/K versus 
X4 for the f.c.c. kattice with A=H/Je=es 10g" 
0.0 and -1.0 from Padé approximants to 
BY(kontamaaan Tne dashed curve is the 


X=] ine for n-s20 265) = 
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For the series in. fugacity we assume the asymptotic 


Form for Y near Zs given by 


NUK hoz) 2 Ceany iT =| 2/2 (kena | OSS Gama 


and analyse the series using Padé approximants to dinY (Kn 5z)/¢4z 
and PaKinz) 7°74. The=second order transition line: in the «tz 
Peonepror the f.c.c. lattice in zero field “iS shown in Figure 
Beoee ane eStimates of y from the analysis are “mot consistent 
except for near the region of what we believe to be the tri- 
critical point. The poles nearest to the origin of Padé 
epproxamants to dinY(K,n,z)/dz are located at Zo and ae on 
weaneor> “Oo but small.° The effect of the pole near -Z, was 
successfully eliminated by transforming the expansion variable 


z to z via the equation 


Desai | (ides Zia) (6.4.4) 


The effect of the transformation (6.4.4) is to leave the 

point Zo invarient and remove the point -Z, to -~ . Estimates 
of Zz. were obtained from the poles of Padé approximants to 
[y(K,h,z]°/" for each temperature considered. Ratios of the 
coefficients B, (Kh) of the transformed series Y(K,h,z) are 
remarkably linear when plotted against 1/2. “ES tamaces sor, 

y(t) from Padé approximants to din¥(Kshozj/dz are given in 


manlerve Ge for the f.c.c. lattice in zero field for selected 
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Figure 6.3 Critical temperature t = Re lzaey 7K versus Zz 


Bor the f.c.c. lattice with n=0.0 from tne 
poles of Padé approximants to 
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Table 6.6 Estimates of y(K) from [N/D] Padé approximants 
to the transformed series dinY(K,h,z)/dz on the 


fuC,C. Maberce itor n=U. 0: I 


[N/D] t=0.,60 CSG. 55 w= 0-90 £=0745 t=0.40 
[1/4] 1.384 isco 0.765 14193 ease 
td a 1.380 1.400 be4i2 Bae. H2 1eah9 
hfe | Vesis 1.401 1.886 ake iahos 
he34 1.433 aoe 1.244 ba2ai00 Teli 6 
2/2) Leo ol deod2 Laas lea a00 Lehkas 


Es 


Baltvessonet <°O0.6°5 “For?’t*> 056 the series sare not well 
benaved and a suitable transformation could not be found to 
smooth the ratio plots. However the estimates of 7Y(K,h) 
from the series in ne snould be the same as from the series 
in Z, since conversion from one series to the otner involves 
only a transformation of variable. Near the tricritical 
point though, series for Y(K,h,z) are expected to behave 
differently tnan the series for Y(K,hyn). 

In the next section we use the series Y(K,h,z) to 


ERy ereuVOcate thei tritriticalspoint inthe itz<p lane. 


Soothe Tricritical Point 

ANT YVcritical point, asi detinedapy Grint itis Clo nye 
iS characterised by the existence of two competing order 
parameters which simultaneously become critical. 

Tne phase diagram for He?-He- mixtures, in the tz 
pene. where, Zz 1s the fugacity of onevof, tne species. (1s Ja 
single continuous line, which nas a tricritical point 
at tL. = Z. andes ie : EO? = Pe the transition is 
from the superfluid to the normal-fluid state of He. When T 
is less. than ibys the transition 1s first order and the dine 
defines the phase boundary of the Het-rich and Hee-rich 
phases. In general all critical exponents are expected to 


have values at the tricritical point different from their 


values at the second order phase boundary. According to the 
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tricritical scaling theory developed by Riedel (1972, 1974) 

and Griffiths (1973) there are two different sets of tricritical 
exponents. Which of these sets is applicable depends on the 
path of approach to the tricritical point and on what variables 
the quantity of interest is expressed in. For our immediate 
purpose st. 1S ‘sufticient..to note that. for fixed kK vand fn it uwe 
eepgoact tne tricritical point using low fugacity series. tne 
exponents are denoted Ops Yeo etc. and have the mean field 


values a, = 1/2 and co - 1. If on the other hand we use the 


t 
low density series the exponents are symbolised bya, vars 


iT) 
ine) 
° 


etc. and have the mean field values oie -1 and ier 


3 Het mixtures by Goellner et al. (1973) 


Experiments on He 
Suggest that tne tricritical exponents have the classical 
values given above. 

The location of the tricritical point using series 
expansions nas been attempted for many different models. 
Systems for which are available expansions for the free energy 
about both the completely ordered and disordered states are 
relatively easily analysed. In such cases, for instance the 
Blume-Capel model as studied by Saul et al. (1974), the first 
order phase boundary may be located by equating the two 
branches of the free energy. Unfortunately, for the spin-z 


XY model the nature of the ground state is entirely unknown, 


so that only expansions about the disordered state may be 


R\y vir i ab 
Se a 
eae | 3 <8 4 
4 Pe pee 
ot hy ghee ote Se oe | 
CW eitae yee ees cia Lae eg ete i Wo Fae + oe a 
he vy ere ne 4 if Tenth Rtas ee : ™ > ae ' 
‘FU AmOLd., ener tail my ib ie been 624 338 ong tl Vee 
i 54 Bi P — hi ty ®) Ma ral 4 rf 
i wai Re , 2 Mie . 
4 ; ¥ WW ' al De? 4 WJ us 
| « “ 4- 
I bed Os a ‘ ee 
7 
fd | of : a pats : ‘ 7 4 
i i i ae 7 P e+ 
r t ¥ i He.ce a ; a ners ‘ 4 ae 
} i Ty ; fe 7 
ah, ; + . me 
rt | ) y" y teak 
+e. 8 has) Vase 13 OF ya zer prey Vom a ns a 
heres Hoe OH Gy 4 eneno 
is _ # her ae era 
<= in q J Al os 
se t ty teh a | re >= * oe py Gating | : bk ed 
obi J0n54 v. winner wel, DI" n't da aa 98 
i 1 * A eG 2 ; , U 
LATIN 4 fz WHY, eng Fee $0n8 @ + det Vexg ote mney, 
Mi zi ; : i) : ere 
ie AAI EIR PIED INS & >t ‘bs sr ebrg rasa par" 
9p sortagent ier a due al Ret 
ize iis ar cite Ps 36 Suet Ye bel oute: | 
) + i’ 
ONE? ae a “Ee iggot 4 dem 
ye oe i * 
eR Te te rau Lieae aaet® at ROLL 
i i re oat auo ad3-4 
/ivuielelns tay : ria eieis siterdil ahi @ g 
. a: a q . aie i i pe a 
- is 
7 7 


ia 


obtained. As yet no precise method has been devised for 
finding a tricritical point using series expansions of 
functions derived from only one branch of the free energy. 

For the spin-z XY annealed site model we nave used 
two methods to locate the tricritical point, botn of which 
assume classical tricritical exponents. The first approach 
we have used is to locate the point in the tz plane where 
eo tee EStimates of Vik), for thé f.c.c. Tattice for 
different values of t, from residues of Padé approximants to 
dinY(K,h=0,z)/dz are displayed in Table 6.7. As can be seen 
from Tables 6.4 and 6.5 the value of y appears to approach 
the mean field value (y = 1) very slowly as a function of 
temperature. Between 0.32 < t < 0.29 y varies only slightly 
from the mean field value. If we assume a tricritical 


temperature t, = 0.31 + .02 then we find using the transition 


t 
line of figure 6.3 that 2, = 0.11 + .02 and from Figure 6,2 
oe O8O<+—~02 

In addition to the fluctuation in order parameter 


the quantity 
w(K,h,z) = zon (K,h,z)/oz (bea) 
is expected to diverge strongly at the tnt Cri tie dlespo nus 


The quantity w(K,h,z) can reasonably be expected to behave as 


a concentration susceptibility because n(K,h,z) is a good 
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Estimates of y(K) from [N/D] Padé approximants to 
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order parameter for the first order transition in the tz 
plane. w(K,h,z) is readily constructed from In HCKenez) as 
a series in z using the relation (6.2.1). We write then the 
esempeotic’ form for wikehsz)oas 775 Z, from below as 


w(K,h,z) ~ ACK,A)ET - 2/z,(K,h)] t (6.512) 


Where we use the exponent notation developed from tricritical 
scaling theory. 
Residues of Padé approximants to dinw(k,h=0,z)/dz 


BOmernte fac.c, lattice at. the point t, = 0.31 yield an 


sR 
mestimate of A, = 0.32 + .10 . This estimate of A, 18 some- 
what lower than the classical value de = 1/2, but because of 


the uncertainty in tne estimate of ty and tne limited lengtn 
of tne series we cannot justifiably conclude that dy jan Ware 
Values of the exponent Ol for the specific heat could not be 
consistentiy estimated. Using the above method Rogiers, 
Dekeyser and Quisthoudt (1975) concluded that x, & 0.8 and 
that the exponents Yeo hy and a, are apparently classical. 
The second metnod we have used involves locating 
the transition line from Padé approximants to [w(K,h,z)] 
assuming Aes 1/2: «a» Tf the curve: obtained nmabhe GtZ plane 
from this method meets the transition line obtained from Padé 
approximants to cv¥(Kh,z)]1/Y with y = 4/3, then we identify 


the point of intersection with the tricritical point. This 
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method was used by Plischke and Betts (1975) to locate the 
tricritical point of the Cheng-Schick model (Cheng and 
Schick (1973)). We have found however, that the line of 
critical points predicted from Padé approximants to 
[w(K,h=0,z)]¢ FOR the fic. Cu lattice fort) <div fens oan 
marginally from the second order line predicted from Padé 
approximants to [Y(K,h=0,z)]°/* displaying in figuresomsr 
For 0.28 < t < 0.31, the poles of the central Padé approxi- 


rd 3/4 


mants to wand Y are the same to witnin 0.03%. 


Consequently nave been unable to make an estimate of ty 
using this method. However, we note tnat the region of 
approximate degenerary of the two critical lines 


(0.28 < t < 0.31) is where tne exponent y 1 


-~ 
t —“~ 


6.6 Summary and Discussion 

The annealed site spin-z XY model nas been studied 
On three cubic lattices. Series expansions up to seventh 
order in magnetic site fugacity and magnetic ion density have 
been derived for the fluctuation in order parameter. Tiese 
series have been analysed in two types of variable, one 


applicable to magnetic systems and the other suited to 


He?-He* 


located in the tz plane. The A-line obtained for n=0.65 in 


mixtures. The tricritical point has been approximately 


the tx plane lies between the A-lines presented by Rogiers, 


Dekeyser and Quisthoudt (1975) for n=0.5 and n=0.8 for the 
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Takagi model. Also we have been able to compare directly the 
second order transition curve for the f.c.c. lattice with H=0 
in the tz plane with that for the Takagi model obtained by 
Rogiers (1974). These two lines agree within the confidence 
limits. 

The estimates of y we have ovtained in the pure 
limit, t=1, do not agree with estimates of y assessed from 
high temperature series for the pure spin-z KY model. “Our 
estimate of y(t=1) is apparently the same as previously 
obtained for the XY model with quenched site impurities. 

As Rogiers, Dekeyser and Quisthoudt (1975) noted, 
the Takagi model (and hence the spin-y XY annealed site model) 
does not predict quantatively the phase diagram of He eHer 
mixtures in the tx plane. The ability of the Cheng=Schick 
model of He?-He® mixtures to match the experimental phase 
diagram (Goellner et al. (1973)) was established by Plischke 
and Betts (1975). Tne Cheng-Schick model allows kinetic 
mobility of tne He? particles sand contains thevconrectfermi 


statistics for these particles. Both of these properties are 


S 


lacking in the Takagi model where the He~ atoms are approximated 


by classical particles with infinite mass. 
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CHAPTER VII 
THE TWO DIMENSIONAL MODELS 


Jeu enicroduction 

Because of the controversy concerning the existence 
of pnase transitions in two dimensional quantum systems, 
we consider the two dimensional quenched and annealed site 
XY models in a separate chapter. 

In Section 2 the exact result of Mermin and Wagner 
(1966), concerning the existence of long range order in two 
dimensional models is briefly explained. Also in Section 2 
tne possible existence of a Stanley-Kaplan type transition 
in sucn models is examined. In Section 3 we present tne 
vesuiasvot Betts, Etliot’ and’ Ditz1anz{(1970) of a high 
temperature series analysis for the XY model on tne 
triangular lattice. Our results from low density series 
for tne quenched and annealed site XY models on the triangu- 
lar lattice are presented in Section 4. A brief discussion 


and summary are contained in Section 5. 


Jeew Exact Results and Conjectures 

Stanley and Kaplan (1966) analysed the nigh 
temperature susceptibility series of Rusnbrooke and Wood 
(1958) for the spin S=1/2, 3/2 and 5/2 two dimensional 


isotropic Heisenberg models. They concluded tnat there 
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was strong evidence to support the existence of a phase 
transition for tne S=3/2 and S=5/2 models but for S=1/2 
the series were not well enough behaved for reliable 
analysis. 

Very ie afterwards nowever Mermin and Wagner 
(1966) proved rigorously for a Varleby.of wiodels ,includind 
the general spin isotropic Heisenberg and XY models, that 
there can be no ferromagnetic or antiferromagnetic order for 
any non-zero temperature. Stanley and Kaplan (1966) then 
made the observation that snould the pair correlation 
functions be only gradually attenuated, then the suscepti- 
bility may become infinite without the onset of long range 
order. Mathematically this means that the susceptibility 
constructed from the sum (2.2.6) becomes infinite without 
the correlation length becoming infinite. Mubayi and Lange 
(1969) have used Green's function techniques with a judicious 
Enoce ort decoupling Seneme to find results for tne s=1/ ic 
isotropic Heisenberg model consistent with the Stanley- 
Kaplan conjecture. Yamaji and Kondo (1973a) using a 
different decoupling scheme arrive at a totally different 
result and soe tainh that there is no Stanley-Kaplan 
Crdysad VON. 

Moore (1969) has derived and analysed series 
expansions for the classical XY and Heisenberg models on 


two dimensional lattices and concludes that there is good 
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evidence for a phase transition. On the other hand Yamaji 
and Kondo (1973b) find, from extended series for the spin-5 
isotropic Heisenberg model, that there is no apparent 
divergence of the susceptibility series. 

Kosterlitz and Thouless (1973) have found a 
possible type of ordered state for some of the two 
dimensional models including the spin ~ XY model. 

Camp and Van Dyke (1975) derived and analysed series 
expansions for the spin ~ XY and Heisenberg models. Their 
results seemingly support tne conjecture of Kosterlitz and 


Thouless (1973) that the susceptibility behaves as 


eect jaeas t > 0 


7.3 Tne Two Dimensional XY Model 

In this section we review some numerical results 
for tne two dimensional spin-5 XY model obtained by Betts, 
Elliot and Ditzian (1970). They derived the high temperature 
series expansion for the Fluctuation in order parameter, 
Vek). for the triangular and square lattices and analysed 
these series using standard techniques and tne method of 
conformal transformations. 

For the series on the triangular lattice, the 
main poles of Padé approximants to dinY(K)/dK are at 
approximately 0.5 + 0.1 i and + 0.6 i. However after 


transforming their series via the quadratic transformation 
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kK = K/(1 - BR) eae) 


the poles closest to the origin are located on the real K 
axis. The effect of the transformation (7.3L 1s bommove 
une points K = + 1/2 uniquely to the points + 1 respectively. 
The series Y(K) obtained after transformation is much better 
behaved than the original series and analysis yields the 
estimates K. = 0.501 + 0.005, with y = 1.50 + 0.02. The 
corresponding value for SN, for tne untransformed series is 

Ke pacror +)0.02. Hence Betts. et “al. conclude wtnatetnesX\ 


model on the trianguiar lattice appears to have a Stanley- 


Kaplan type of phase transition. 


7.4 Analysis of the Density Series 

In this section we present tne results of the 
analysis of series for the fluctuation in order parameter 
eneche triangular lattice in terms» of the magnetic site 
density for the quenched and annealed site XY models. Since 
we are concerned witn obtaining evidence for or against the 
existence of a phase transition we will examine in detail 
one series only. The series we use is tne fluctuation in 
order parameter as a series in density p, for tne quenched 
Site case at K = 0.65 ~ Ko . A preliminary investigation of 
the quenched site series for various values of K indicated 


by extrapolation that at K = 0.65, Po would be approximately 
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1.0 . Scaling theory predicts that the critical behavior of 
this series should be the same as for the nigh temperature 


series for thevpure XY model. The ‘series we Study below is 


Z 3 


¥(p,k=0565) = pu+ 1,8841255.-.po401.93644915. Ip 
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Tne nigh temperature series for the free energy 
f and fluctuation in order parameter Y for the quenched 


Site XY model on tne triangular lattice are: 


TA PRR }PSs In2e + 4 0* Ke ° 5 p °K? 4 (~55 po “ 2 ps : 3 p4)kt 
“i (-7 p> - 5 pt + 2 pK + (qua P* + co >? x aN af 
rea p> + a p°)K° 
+ GS p> + ae p4 3 fn p? - 2 »° + op’ )K! 
‘: (74.2) 


and 


Wip.K)/p = 1 + 3pk + Fy ph? + (-L p - 24 p? + 17h p23 
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Themseries (7.4.2) and (7.4.3) are of insufficient length 
to analyse with confidence. The series of Betts and 
co-workers referred to in the previous section are several 
terms longer. 

We present now the results of the analysis for 
the low density order parameter fluctuation series (7.4.1) 
Estimates of Po from poles of Padé approximants to 
dinY(p)/dp are given in Table 7.1. The estimates of p. 
presented in the table are not the poles closest to the 
origin however and consequently do not give consistent 


estimates of the location of tne physical singularity. In 


each case the pole closest to the origin is on the negative 


rea leaxicwat about p'@ -0. oo  Ineonder “to moyvercciirs 


singularity outside the radius of convergence of tne series 
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ape 7.) Estimates of p. and, in parentheses, y from Padé 
approximants to the series for dinY(p)/dp. A 
missing entry means that all poles of that 


approximant are unpnysical 


D N=0 N=1 N=2 N=3 
2 059737 0.9348 1.0153 TeOZZ0 
(0.688) (1.076) (i370) (1.406) 
3 0.9945 = WAO22u, 
(2324) (445) 
4 1.1840 1.0579 
(4.756) (1.699) 


| 
Cote —~ ob andes © 
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we make the transformation 


Dea b/ (2 > p/0. 9) (P7pe Area) 


The effect of the above transformation is to leave the 
point p=0.8 invarient and to remove the DOINtsp = =U. catc 
pr-o. 

Table 7.2 shows, for Padé approximants to 
dinY(p)/dp, the closest poles to the origin and their 
residues. The low order Padé approximants are included 
for completeness only, but it is interesting to observe the 
consistency of estimates of Po throughout the table. We 
conclude that p. = 0.885 + 0.005 (p=0.990 + 0.005) and 
Ves s23 + 0.05 by taking an average from the higher degree 
central Padé approximants. Our estimate of y. differs 
considerably from tne estimate of y found by Betts et al. 
(1970). Ff indeed y#y, then a violation of universality 
icmampliaed, just as for the three dimensionalecase. sWe can 
however analyse series for alee fOr Varvous -choices. of 
Y, in order to determine an optimum value for y from the 
consistency of the poles of the Padé approximants. 
infortunatety: we find that the canalysis7otssenives =a or 
v(p) '/Y are insensitive to the choice of y. In Table 7.3 
the poles of Padé approximants to Y(p) tf? are displayed. 


The central Padé approximants give an average of 


i> , iy, 
Cae fd oan 
¥ | ! ¢ 
i 7 
ce : " 
i y i i 
: , i ly 
Fi = we hy 7 
; ” ~ ny: he) 
‘oe toe uM ve ee 
v 7 hice P an a rs a ] a B , ; 7 - 
, ate, yer Ri , Baa f ‘ ae ae 
a Ac ; ' i : yo ; . ng re) } - {ee 
rome et eee ee Fok i, a BSc... 
¥ ae z f } , (a 
ri ‘ Prive eT PO ers ’ ’ ae Pe eee | 2 
j Cae 7 vo i ee P ’ re OOS ety ee ba Am 
ea Me ED i. Te ar neh ve 
‘ aan an Le: 
A Oe) a 
} Mig, ! 4 
os 3 « 1 men ; ae a oe Oe 4 y i | > fan Ds 
F P * t 
' ‘ 7 
i ro. 
‘ ae en ‘ % » Ph le 
‘ P ‘ 5 cay... 
\ a oe ri 7 
' ‘ait 5 Wiest a ae ae m4 j r >, 2 
4 ra - i i 
Ny : 
7h -— i \ wo 
> ‘ 4 Lue ckin® 
a re ey { b 
: ‘es YW OR Se f 2% 
r d if ;| 
I 
bie 
rm a oe > + » ie = 


, ‘ ( of 
" ror 
’ f o k d 1 - r 
(‘ts i hie aE: Ota WG) na 
f in uf 
' :; . 4 
Y ; } , 2 ‘ >) 
7 bt  * ih i ¥ | tb) Cy 
hy bi ow : Bie rer k . ’ * iy ys 


¢ ] ; r a a | P vr 22) ro h : ee | ot 
et bie ae a ae Tok yis:coe a ee 
i i Fe ty ike , A ) rl 


rovepeeeeta se nab) ed nom xoiqds She 


.s 


+ esis 


we a 
Prd 
f 


3 A Mita yy ov estman inicbeas 


- 


OES ae oa *, 
by Mino 1b GH ee lait 5. nat aa ‘oi 
an a Beg oan y 
ped % Oo <8 [ ea . ans a he Q? 
aes 3 i : i 5 Ns 
= any yi a Hae’ ‘Bark wy 
me yt 


oe . 7 


= 


y. end oa 


i, sh on £x0% 


i; i 
" . 5 TiN, 1 ae 


rZ26 


Babiecwd-c Estimates of ee and, in parentheses, y from [N/D] 
Padé approximants to the transformed series 


dinY(p)/dp 


D N=0 N=] Wed N=3 N=4 


] 02.95.03 ONS 0. 690.1 0.8842 0.8845 
(1.48) ea) (1.43) (4 0p) (1.24) 
2 0.9151 | 0.8583 0.8816 0.8845 
Rae 2775) (1.04) Cie 2y a2 by) 
3 8.5935 0.8818 0.8868 
fern) (eteea2)) (iz) 
4 0.886] 0.8849 
reetcl:) (1.24) 
5 0.8850 
eeiiC:) 


{em 6 ee ee 


Table 7.3' Estimates of Po from poles of [N/D] Padé 


approximants to V(p)4/° 


D N=0 N=] N=2 N=3 N=4 


] Oo 17 0.8673 0.8832 0.8860 076552 
2 0.8747 0.8881 0.8865 0.8859 0.8860 
3 0.8864 0.8866 0.885] 0.8856 

4 0.8866 0.8864 a peveken®) 


5 0.8859 0.8856 


Wee 


0.8857 


a8. © 
Dr = 0.8856 + 0.0006, which corresponds to PL 0.993) Oe 
The precision of this estimate for Be is misleading though, 
because poles of Padé approximants to ¥(p)2/9 give almost 
the same estimate for Pes with only slightly less precision. 
Tne series for the fluctuation in order parameter 
Y(n.) in powers of the magnetic site density n,, on the 


Epranguiar lattice: for K = 0.65 1s 


2 3 
Y(n_) Sere) 20202 02.088 en tec 2005 Ae 
4 5 
40273523268". -hoin st. 2. 4860704) Wea ne 
6 ii 
ED CeaOly 72) OO Gees tet nee aE OOO 0/1) meee ree Ne 
Hae) C7e Geos) 


The results of the analysis of the series (7.4.5) 
are very similar to that presented above for the series 
(7.4.1) and need not be presented in detail. 

We were unable to obtain more than a small part 
of the line of transition points in the temperature wdencucy 
Diane.) inthe region 0.9 <-t 1.0.7 wihene  t=h7 010s) aaon 
the quenched site case the critical probability Po is 


approximately given Dy 


po = (t + 12/2 (77 aes) 
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stim ye - 
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sia ot nee. 


Vor 


a similar formula nolds for the annealed site case. . The 


mean field result is 


sOurne. inttial slope for the line of critical point differs 
greatly from tne mean field result. This contrasts with the 
results found for the three dimensional models, where the 


initial slopes were very close to the mean field values. 


7.5 Summary and Discussion 

In the preceding section, evidence in support of 
a Stanley-Kaplan type phase transition for the XY model on 
the triangular lattice has been presented. Our results 
corroborate those of pages. Yoni ts and Ditzian (1970) in 
respect to tne existence of a phase transition and for the 
pure XY model indicate tnat Ko = 0:67 (4°0.04." On the other 
hand the estimates of y we obtained do not agree with the 
Valge y= 1.50 found by Betts et al. 

Recently, application of ‘real space renormalisation 
group theory to the pure XY model has yielded confusing 
results. These results are reported briefly below. Rogiers 
and Dekeyser (1976) find, using a three spin cell on the 
triangular lattice, that there is a transition for the XY 
model at about K.~ 0.8 . Betts and Plischke (1976) find 


however, using a five spin cell on the square lattice that 
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tmereris no fixed pointiliees mo iplase transition): tor 
Wee... ALSO Rogiers and Betts (1976); usingsa seven spin 
eetl-on tne triangular lattice, find Sac Lot) =e Rn emmesults 
of both Rogiers and Dekeyser (1976) and Rogiers and Betts 
(1976) are suspect because negative values are found for the 
exponent 6, which is thermodynamically impossible. 

In short the evidence for a Staniey-Kaplan 


transition for the two dimensional XY model is inconclusive. 
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CHAPTER VITI 
SUMMARY OF RESULTS AND DISCUSSION 


The quenched and annealed site diluted spin-z XY 
models have been studied in some two and three dimensional 
Bravais lattices. We have used the finite cluster theorem 
(outlined in Chapter IV) to construct expansions about the 
magnetically disordered state for the free energy and the 
fluctuation in order parameter. 

For the quenched site model, expansions have been 
constructed up to degree seven in both magnetic site 
concentration, p (with De Rmeae ane as a parameter) and 
in inverse temperature (with p as a parameter). The natural 
expansion variable for the annealed site model is the 
magnetic site fugacity z. Expansions up to degree seven in 
z have been derived and these series have also been 
converted to series in magnetic ion concentration, nes 

The results of the analysis of the series expansions 
for the quenched site model for the three dimensional 
Lattices are given in Chapter V.. In brief, we have found 
for each ieenices thésline “of critical points” tn. the 
temperature-density plane fairly precisely. However, the 


apparent behavior of y(t), the susceptibility exponent, 


(where t = T(1)/T(p)) has some disconcerting properties. 
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Firstly y(t=1) = 1.22 + 0.02 does not agree with the 

value y(p=1) = 1.333 + 0.002 obtained from the nigh 
temperature series approach. This result contradicts the 
scaling hypothesis, which asserts inter alia that 

y(t=1) = y(p=1). Also we have found that Y(t) apparently 
varies with reduced temperature t. This contradicts the 
universality principle which requires that y(t) be indepen- 
dent of t. Very similar results have been found for the 
annealed site caSe. 

For the quenched site case it can be argued’ that 
the concentration p is not a thermodynamic variable and 
hence the phenomenalogical theores of scaling and universality 
need not apply. This argument is weakened by the fact that 
the annealed site model shows the same behavior for y(t), 


Since the concentration nia for the annealed site model is a 


— 


thermodynamic variable. 


The phase transition curves for the annealed site 
model are given in Chapter VI. The interpretation of the 
annealed site model as a model of He>-He? mixtures is also 
given in Chapter VI along witn out attempts to locate the 
tricritical point. Tne model does not appear to give a good 


, 3 4 
quantitative description of tne phase diagram for He -He 


mixtures. 
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Tne two dimensional diluted XY models are examined 
in Chapter VII and our results seem to indicate the presence 
of a Stanley-Kaplan type phase transition. 

Possibly one of the more interesting results we 
have found is that the nature of the phase transitions for 
the three dimensional lattices is uncomplicated by the 
introduction of impurities. However if there were, for 
instance, no precise critical density for a given temperature, 
but instead a small interval of critical densities, then 
analysis of expansions about the disordered state only might 
be unable to reveal the distribution of singularities. 

The only important issue we have left unresolyed is 
the discrepancy between our deduced behavior of y and the 
predictions of scaling and universality. This «dtfficul ty 
cannot, in the foreseeable future, be resolved using series 
expansions because of the size of the eigenvalue problems 
involved in extending our series. Possibly a more useful 
approach would be to examine the problem from a general 
thermodynamic viewpoint. 

The weak universality hypothesis could be tested by 
constructing series expansions for pair correlation fem.ctions:. 


The ratio of the exponents y/v should be invarient along the 


line of critical points. 
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APPENDIX A 
EXPECTATION VALUES FOR FINITE CLUSTERS 


In order to utilise the finite cluster expansion as 
proposed in Section 4.5 we need the partition function and 
expectation, values: for teach -finitbe cluster. 

If we already have the Hamiltonian for a cluster 
eyeot Say. V vertices, in matrix form, then it sis sa Simple 
foecers tO ind the partition ufuncti on of.€. Al Lithat sis 
required is to diagonalise the Hamiltonian to find the 
energy spectrum eri ya? naa een (The XY spin-z mode] 
Mmas.a2 States per site and hence for-a v. vertex cluster 
there are 2” states in all.) The Dameit ton sfunct..on) is 


tnen 


O(C) ="Tr expC-BH) = exp(-Be,) (CAeain) 


ros) pe) 
ibaa 
— 


where £8 = 1/kpT . The sexpectation value (C) forthe 
V 


operator y= 2 Se Si for the cluster C is given by 
i, j=l 
y = Tr y exp(-6H) 
oV 
ey ae LY, exp(-Be, ) GANS 2) 
Vas a 


where yo, 1s the &£th diagonal element of the matrix for y 


in the representation where the Hamiltonian is diagonal. 
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The Hamiltonian 


pee Sy (o%o% + oo%) - mi x s? (A.3) 
<j J J j 1 


then must be constructed for each finite cluster. The most 
Straight forward approach is to write out the Hamiltonian 
in the direct product representation, where for a cluster of 
v spins IA is }qivenaby Glut = :ta4-) wor 1. peWhere whew fac 
tiesa component of the’ Pauli spin» vector and is in ther ith 
position in the sequence of direct products of the 2x2 
identity matrix I. 

By constructing the matrix for x $3 in the direct 
preduct representation it can fred vaeeay be seen how the 


xv block diagonalises into matrices 


interaction part of H 
corresponding to the same a2 eigenvalue. This approach is 
best iliustrated by example. 


Consider the labelled cluster 


2 3 


y 


for which we have 2 S4 = diag(4,2,2,0,2,0,0,-2,2,0,0,-2,0, 
j 
-2,7-2,-4) in the direct product representation, where since 


zi sf has non-zero elements on the diagonal only, we write 
cera 

7 ; 
down only the diagonal elements. It is obvious then that if 


we exchange rows and columns 4 with 9 and 8 With 13. the 
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é Z ; . 
matrix & 3 will be in blocks of equal eigenvalues. If we 
j 


also exchange tne same labels in the matrix representation 
INT 


x 
Oh al = x (aro ovoy), the result is the matrix 
Sh mane a a 


acid fo 1 0700. O@ieniral 
fon Grae © / 10 0.1.0 0 \°@ {1-000 saeco 
t= 180.070 10 7OmOre | Gor 4 0, 0d 
Neti ie! Ql SOKO O Ny 12080807 
\ O,.G te Op Olly 
Svante. Lop lee pated! 


Thus the problem size is reduced from a 16x16 matrix to that 
of a 4x4 and a 6x6 matrix. These matrices are readily 

diagonalised to give the set of eigenvalues 0, 0, 0, +73, +1, 
fete 0 0, +v3, 0 and in zero field ithe partition epunaeuen 


for this cluster is 


Q = 6 + 4Coshv3K + 4CoshK + 2Cosh2K 


where K = BJ. 


Tne above example essentially illustrates the 
computational technique used to construct the partition 
functions of each cluster. To convert the y matrix from 
the direct product representation to the energy representa- 


tion requires that we form the matrix 
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where v is tne matrix of eigenvectors obtained by diagona- 
lising the Hamiltonian. Standard routines were used to 
diagonalise the Hamiltonian and the entire process was 
computerised. 

For symmetric clusters a further block diagonalisa- 
tion can be achieved by writing the Hamiltonian in a 
Symmetrised basis. However, the advantage of tnis refinement 
is offset by the relatively low symmetry of most of the 
graphs and tne fact that the transformation to the 


symmetric basis is unique for each cluster. 
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APPENDIX B 
LATTICE CONSTANTS AND PERIMETER COUNTS 


The following table contains the strong lattice 
constants and perimeter counts for all the connected 
clusters of up to and including six vertices. These data 
are complete for tne face centred cubic ( FCC ), the body 
Gentred cubie (BCG ie) : the simpleecubienGSCu, ) sethe 
triangular ( TRI ) and the simple quadratic ( SOU ) 
lattices. The graphs are in dictionary order and sub- 
headings give the number of vertices and number of bonds. 
Each graph is represented by an adjacency matrix A such 
that the matrix element As j is one if there is a bond 
between the ith and jth sites and zero otherwise. For 


example the graph 
3 4 


A/ 


2 ] 


(number 7 in the list) has adjacency matrix 


FO: site ok 
Nig 8 Aol ball 
} ye ad 
12 020 20; 


147 


a 
ae 
ve . 4 
a nih Buoy sot ombaaie: er 
; my ee sMAP pated ont: naecaa a 2 8 
J : i = oe ie ¥ ‘ in 
5 an | “ie y 
AS A> kD Wyded Weranes ooctheebinel wast 
b= , je ei] thee, ot ARTE BAP e4 te <a ae ahaa 
yy) SD) aig ies wih’ a fateh at orl sbeees joan aa 
¥, ery RY pahee .) oe ty G } j 4 ih Ki 4) 4i9 anew ante ell 
ied HAM AE asatoie) Fo Aes aura wit avto § : 
. . / are ac as | F on 
as a | 472 POS IK fade + > “cd . tra ri9esnaey, cael hi ica 
| Neds 6 eh aera, FF 900 et) yA, Pee: pewvensi ie 
14 -  Seteretia ovat bos..carste. ase cae ae 
, pir ee 


' 


sitiam \ votre tha seit ‘(oath wea € 


’ a . 


7 , 
x we ny ep ind, e. Pin i‘ 5 
i J en rae +) 


148 


Because this matrix is symmetric and all the diagonal 
elements are zero the graph is adequately represented by 
the above diagonal elements collected column by column 


Pomnorm the ‘vector, “A= 1 1 1 7 0 0:6 


An example of how these data are to be read is 
given below. 
Consider the graph with four vertices and four bonds 


for which the following data are listed 


7 2 111100 
FEC: ie: 48(27) 72(28) 
TRI: Tas Nelle) 


This graph is the seventh in the list and has a symmetry 
group of order two (i.e. symmetry number of two). The next 
set of numbers is the adjacency matrix in the form given 
above (in fact this is the graph used in the above example). 
The istrong lattice constant for the seventh graph is 120 on 
the f.c.c. lattice and is embeddable 48 times with perimeter 
27 and 72 times with perimeter 28. On the triangular lattice 
the graph can be embedded 12 times with a perimeter of 11. 
Note that if a graph is not strongly embeddable in a 


particular lattice, then no data for that lattice appear. 
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1 VERTICES } BONDS 
1 BZ 0 
FCC: 1; 1 (12) 
BCC: 1; 1(¢ 8) 
SCU: 1; ao a) 
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2 2 eet 
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